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Abstract 

A  developable  surface  can  be  formed  by  bending  or  rolling  a  planar  surface  without  stretch- 
ing or  tearing;  in  other  words,  it  can  be  developed  or  unrolled  isometrically  onto  a  plane. 
Developable  surfaces  are  widely  used  in  manufacturing  with  materials  that  are  not  amenable 
to  stretching.  A  ship  hull  design  entirely  composed  of  developable  surfaces  would  greatly 
reduce  production  costs  of  that  hull. 

This  thesis  describes  a  new,  user-friendly  method  of  designing  developable  surfaces  with 
a  B-spline  representation.  First,  it  expands  the  work  of  Aumann  in  designing  developable 
B-spline  strip  surfaces  whose  directrices  lie  on  parallel  planes.  The  computer  program 
developed  can  assist  in  designing  surfaces  from  curves  of  any  degree  with  any  number  of 
segments,  and  includes  a  test  for  regularity.  Second,  a  new  method  is  developed  which 
permits  the  design  of  developable  surfaces  with  general  three-dimensional  space  curves  as 
directrices.  The  computer  program  permits  design  of  degree  (3-1)  B-spline  developable 
surfaces  with  two  patches  using  a  minimization  process.  The  basis  is  provided  to  extend 
this  to  surfaces  with  more  patches  and  higher  degrees.  A  test  is  provided  to  ensure  that  the 
minimization  process  results  in  a  developable  or  nearly  developable  surface.  This  method 
is  extended  to  include  special  interesting  cases  such  as  triangular  degenerate  patches  and 
surfaces  with  a  planar  resulting  directrix. 

This  thesis  also  treats  common  differential  geometry  properties  such  as  lines  of  curvature 
and  geodesies  that  are  useful  in  the  design  and  manufacturing  process.  Lines  of  curvature 
are  vital  in  the  forming  process,  since  the  planar  shape  must  be  placed  so  the  rollers  are 
parallel  to  the  In  • ;  of  ~;ero  curvature.  As  an  inflection  line  greatly  affects  the  forming 
process,  a  method  is  aescribed  to  determine  the  inflection  line  in  advance;  several  properties 
of  developable  surfaces  related  to  inflection  lines  are  also  described. 

Straight  lines  on  a  plane  map  to  geodesies  on  a  developable  surface.  This  fact  is  of 
assistance  in  many  aspects  of  the  use  of  developable  surfaces,  including  layout  and  quality 
control.  A  method  is  described  to  calculate  geodesies  on  a  developable  surface  as  an  initial 
value  problem  instead  of  the  more  complicated  boundary  value  problem. 

Engineering  examples  are  provided  for  each  topic,  and  a  small  boat  hull  is  designed 
using  the  methods  described.  In  addition,  recommendations  for  further  research  are  made. 

Thesis  Co-Supervisor:  Nicholas  M.  Patrikalakis,  Professor  of  Ocean  Engineering 

Thesis  Co-Supervisor:  Takashi  Maekawa,  Lecturer  and  Research  Scientist 

Thesis  Reader:  David  C.  Gossard,  Professor  of  Mechanical  Engineering 


Acknowledgements 

I  would  like  to  express  my  sincere  appreciation  for  the  guidance  provided  by  Dr.  Takashi 
Maekawa  and  Prof.  N.  M.  Patrikalakis.  I  would  additionally  like  to  thank  Prof.  D.  C.  Gos- 
sard  for  his  assistance,  Mr.  Stephen  Abrams  and  Mr.  Fred  Baker  for  their  assistance  with 
computer  work,  Dr.  Xiuzi  Ye  for  his  feedback,  and  the  members  of  the  Design  Laboratory 
including  Dr.  Wonjoon  Cho,  Mr.  Todd  Jackson,  Mr.  Guoling  Shen  and  Mr.  Guoxin  Yu  for 
their  support. 

Special  thanks  to  LT  Allan  Andrew,  LT  Kurt  Crake,  LT  Tom  Laverghetta,  LT  Chris 
Levesque,  LT  Casey  Moton,  and  LT  Tom  Trapp  for  their  camaraderie  and  support,  and  to 
LCDR  Mark  Welsh  and  CAPT  Alan  Brown. 

This  work  was  funded  by  the  United  States  Navy. 


Contents 


Abstract  3 

Acknowledgements  4 

Contents  5 

List  of  Figures  7 

List  of  Tables  9 

1  Introduction  11 

1.1  Background  and  Motivation 11 

1.2  Research  Objective 12 

1.3  Thesis  Organization 12 

2  Review  of  Curves  and  Surfaces  13 

2.1  Introduction 13 

2.2  Basic  Theory  of  Curves    13 

2.3  Basic  Theory  of  Surfaces 18 

3  Curve  and  Surface  Representation  25 

3.1  Introduction 25 

3.2  Bezier  Curves  and  Surfaces 25 

3.3  B-Spline  Curves  and  Surfaces 29 

4  Review  of  Differential  Geometry  Properties  of  Developable  Surfaces  37 

4.1  Introduction 37 

4.2  Ruled  and  Developable  Surfaces 37 

5  Design  With  B-Spline  Developable  Surfaces  47 

5.1  Introduction  and  Literature  Review     47 

5.1.1  Developable  Bezier  Strips 47 

5.1.2  Duality  Between  Points  and  Planes 48 

5.1.3  Nonlinear  Representation 48 

5.1.4  Development 48 

5.2  B-Spline  Developable  Surfaces  with  Directrices  in  Parallel  Planes 49 

5.3  Developable  Surfaces  with  3D  Directrices     54 

5.3.1  Developable  (n,  1)  Bezier  Surfaces 54 

5.3.2  B-Spline  Developable  Surfaces     59 


5.3.3  Accuracy  of  Resulting  Developable  Surfaces 62 

5.3.4  Developable  Surfaces  Constrained  by  Planes 63 

5.3.5  Degenerate  Developable  Surfaces 64 

5.4  Development  of  a  Developable  Surface  onto  a  Plane     66 

5.5  Examples 68 

5.5.1  Example  of  a  Developable  B-Spline  Strip  Surface 68 

5.5.2  Example  of  a  Developable  B-Spline  Surface  with  3D  Directrices  ...  69 

6  Geodesies  on  Developable  Surfaces  73 

6.1  Introduction 73 

6.2  Formulation 73 

6.3  Examples 77 

7  Lines  of  Curvature  on  Developable  Surfaces  81 

7.1  Introduction 81 

7.2  Inflection  Line     81 

8  Engineering  Example  89 

9  Conclusions  and  Recommendations  95 

9.1  Conclusions 95 

9.2  Future  Work    95 

Contents  97 


List  of  Figures 


2-1  Tangent  Vector 14 

2-2  Osculating  Plane 15 

2-3  Space  Curve .  17 

2-4  Curvature  Vectors     20 

3-1  A  Cubic  Bezier  Curve  with  Control  Polygon      26 

3-2  The  deCasteljau  Algorithm 27 

3-3  A  Bezier  Surface  with  Control  Net 28 

3-4  Continuity  Conditions 29 

3-5  A  Cubic  B-spline  Curve  Segment  with  Control  Polygon     30 

3-6  Convex  Hull  Property  for  a  Cubic  B-Spline  Curve     31 

3-7  The  deBoor  Algorithm 32 

3-8  Knot  Removal 34 

4-1  A  Ruled  Surface    38 

4-2  A  Developable  Surface  with  Tangent  Plane  along  a  Ruling 38 

4-3  An  Envelope  of  a  Family  of  Curves 40 

4-4  Edge  of  Regression 41 

4-5  A  Surface  with  Folds  and  a  Developable  Surface     44 

5-1  A  Developable  Surface 48 

5-2  Design  of  Developable  B-Spline  Strip  Surface 50 

5-3  Effect  of  Simple  Bounds  on  Developable  Surfaces 57 

5-4  Design  of  Developable  B-Spline  Surfaces 61 

5-5  Planar  Constraint  on  Developable  B-Spline  Surface  Design 64 

5-6  Triangular  Degenerate  Cubic  Bezier  Patch 66 

5-7  Development  of  B-Spline  Developable  Surface 69 

5-8  Automobile  Windshield,  a  Developable  Degree  (4-1)  B-spline  Strip  Surface  70 

5-9  Ship  Stack,  A  Developable  Degree  (3,1)  B-Spline  Surface  with  \K\  <  2.19  x 

10~7 71 

6-1  Geodesies  on  a  Closed  Surface 74 

6-2  Geodesic  on  a  Degree  (3,1)  Developable  Surface 75 

6-3  Geodesic  on  a  Degree  (4,1)  Developable  Surface 78 

6-4  Geodesic  on  a  Degree  (4,1)  Handkerchief-like  Developable  Surface 79 

7-1  Developable  surface  and  its  control  polyhedron  with  inflection  line 86 

7-2  (a)  Lines  of  curvature  of  developable  surface  with  inflection,  (b)  Magnifica- 
tion near  inflection  line 86 


7-3     (a)  Lines  of  curvature  on  perturbed  surface  £  =  0.02.  (b)  Magnification  near 

u=0.57 87 

7-4     (a)  Lines  of  curvature  on  perturbed  surface  £  =  0.08  (b)  Magnification  near 

u=0.57 87 

8-1     Planing  Boat  Terminology 89 

8-2     Boat  Composed  of  Developable  Degree  (3-1)  B-Spline  Surfaces 91 

8-3    Boat  Bottom  View 91 

8-4     Boat  Surfaces  Developed  onto  a  Plane 92 

8-5     Geodesies  on  the  Forward  Side  Section 93 

8-6     Lines  of  Curvature  on  Boat  Surfaces  (Solid  Lines  are  Maximum  Curvature 

Lines,  Dashed  Lines  are  Minimum  Curvature  Lines)     94 


List  of  Tables 


6.1  Control  Points  for  a  Degree  (3,1)  Developable  Surface 77 

6.2  Control  Points  for  a  Degree  (4,1)  Developable  Surface 78 

6.3  Control  Points  for  a  Degree  (4,1)  Handkerchief-like  Developable  Surface    .  .  79 

7.1  Control  Points  for  a  Degree  (3,1)  Developable  Surface 85 

8.1  Target  Curve  Data 90 

8.2  Developable  Surface  Data,  Starboard  Side 90 


10 


Chapter  1 


Introduction 


1.1     Background  and  Motivation 

A  ruled  surface  is  a  curved  surface  which  can  be  generated  by  the  continuous  motion  of 
a  straight  line  in  space  along  a  space  curve  called  a  directrix.  This  straight  line  is  called 
a  generator,  or  ruling,  of  the  surface.  Developable  surfaces  are  a  subset  of  ruled  surfaces 
which  have  the  same  tangent  plane  at  all  points  along  the  generator.  A  developable  surface 
can  be  formed  by  bending  or  rolling  a  planar  surface  without  stretching  or  tearing;  in  other 
words,  it  can  be  developed  or  unrolled  isometrically  onto  a  plane.  Developable  surfaces  are 
also  known  as  singly  curved  surfaces,  since  one  of  their  principal  curvatures  is  zero. 

Developable  surfaces  are  widely  used  in  manufacturing  with  materials  that  are  not 
amenable  to  stretching.  Applications  include  the  formation  of  ship  hulls,  ducts,  shoes, 
clothing  and  automobile  parts  such  as  upholstery  and  body  panels  [14]. 

In  shipbuilding,  developable  surfaces  are  shaped  using  only  rollers  or  presses.  Heat 
treatment  is  then  used  only  to  remove  distortion  induced  by  welding  or  other  means.  Dou- 
bly curved  surfaces,  on  the  other  hand,  must  be  heat  treated  after  rolling  to  induce  the 
additional  curvature.  The  heat  treatment  is  normally  done  by  hand,  by  a  skilled  artisan 
with  years  of  training  to  achieve  the  correct  amount  of  bending.  This  is  an  extremely  time 
consuming,  labor  intensive  and  thus  expensive  process. 

According  to  Avondale/IHI  Shipbuilding  Technology  Transfer  data  for  a  tanker,  only 
15.1%  of  the  curved  plates  in  a  ship  hull  are  singly  curved,  while  65.8%  of  the  plates  are 
doubly  curved,  requiring  roller  and  heat  treatment  processes  [28].  An  intensive  effort  to 
increase  the  amount  of  developable  surfaces  in  the  hulls  of  merchant  ships  at  Burmeister 
&  Wain  Shipyard  [34]  has  resulted  in  a  reported  20%  reduction  in  manhours  required  to 
produce  a  hull.  Designing  a  ship  entirely  of  singly  curved,  or  developable,  surfaces  would 
reduce  construction  costs  even  more. 

Recently  researchers  in  Computer  Aided  Geometric  Design  have  been  quite  active  in 
investigating  the  representation  of  developable  surfaces  in  terms  of  Non-Uniform  Rational 
B-Spline  (NURBS)  surfaces  or  a  special  case  of  NURBS  called  Bezier  surfaces  [1,  3,  14, 
29].  NURBS  curves  and  surface  patches  are  the  most  popular  representation  method  in 
CAD/CAM  due  to  their  generality,  excellent  properties  and  incorporation  in  international 
standards  such  as  IGES  (Initial  Graphics  Exchange  Specification)  and  STEP  (Standard  for 
the  Exchange  of  Product  Model  Data).  Thus,  it  would  be  beneficial  to  design  developable 
surfaces  using  a  B-spline  representation. 
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1.2  Research  Objective 

Recently,  efforts  have  been  underway  to  revitalize  commercial  shipbuilding  in  the  United 
States.  Several  professors  at  MIT  are  involved  with  developing  cost  saving  methods  in  ship 
fabrication  to  ensure  that  shipbuilding  in  the  US  will  be  competitive  in  the  world  market. 
This  research  is  in  support  of  one  of  those  efforts. 

The  main  goal  of  this  research  is  to  develop  a  user-friendly  method  of  designing  devel- 
opable surfaces  with  a  B-spline  representation.  The  effort  is  then  extended  to  address  some 
common  differential  geometry  properties  that  will  be  useful  in  the  design  and  manufacturing 
process. 

The  ultimate  goal  is  to  provide  a  method  to  design  a  complete  ship  hull  from  developable 
surfaces  and  to  generate  cutting  and  bending  information  in  a  format  that  is  user  friendly 
for  both  the  engineer  and  the  worker.  Although  this  thesis  does  not  go  that  far,  it  takes  a 
major  step  toward  this  goal. 

1.3  Thesis  Organization 

In  this  thesis,  Chapters  2  through  4  review  basic  differential  geometry  properties  and  intro- 
duce the  concepts  of  developable  surfaces.  They  also  briefly  review  the  Bezier  and  B-spline 
representations  of  curves  and  surfaces.  The  information  introduced  in  these  chapters  will 
be  used  throughout  the  thesis. 

Chapter  5  describes  a  new  user  friendly  method  for  design  of  developable  surfaces  in 
B-spline  representation.  Section  5.1  reviews  the  current  literature  on  developable  surfaces. 
Section  5.2  describes  the  design  of  strip  surfaces  that  are  constrained  between  two  parallel 
planes.  Section  5.3  describes  the  design  of  surfaces  with  directrices  that  are  space  curves 
and  treats  special  cases  such  as  triangular  degenerate  patches  and  surfaces  with  a  planar 
resultant  curve.  Section  5.4  describes  the  unrolling,  or  development,  of  the  surfaces  into  a 
plane.  Section  5.5  gives  some  engineering  examples. 

In  Chapter  6,  a  geodesic  on  a  developable  surface  is  found  as  the  solution  to  an  initial 
value  problem  rather  than  a  two  point  boundary  value  problem.  The  geodesies  can  be  used 
in  ship  design  for  laying  out  butts  and  strakes,  among  other  uses.  In  Chapter  7,  the  lines  of 
curvature  on  a  developable  surface  are  analyzed  and  the  line  of  inflection  is  defined.  These 
concepts  are  required  for  the  bending  of  steel  plates  into  developable  surfaces;  the  steel 
must  enter  the  rolls  in  a  direction  parallel  to  the  lines  of  zero  curvature,  and  cannot  be 
rolled  past  a  line  of  inflection  [33]. 

A  small  boat  was  designed  using  the  methods  described  in  this  thesis.  The  boat  is 
presented  in  Chapter  8  as  a  practical  engineering  example. 

The  thesis  concludes  with  Chapter  9  which  includes  recommendations  for  further  re- 
search. 
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Chapter  2 

Review  of  Curves  and  Surfaces 


2.1     Introduction 

This  chapter  reviews  the  basic  theory  of  curves  and  surfaces  including  such  topics  as  the 
Serret-Frenet  formulae,  the  first  and  second  fundamental  forms,  curvatures  and  geodesic 
curves  which  will  be  used  in  derivations  of  many  formulae  in  the  later  chapters.  The 
information  in  this  section  follows  the  derivations  in  classical  texts  on  differential  geometry 
and  geometric  modeling  such  as  those  by  Kreyszig  [26],  Struik  [41]  and  Faux  and  Pratt  [11]. 


2.2     Basic  Theory  of  Curves 

Throughout  this  thesis,  curves  are  represented  in  a  parametric  manner,  r(u),  as  a  function 
of  one  parameter,  u,  that  lies  within  a  closed  interval  u\  <u  <  U2-  The  curve  is  a  mapping 
from  a  one  parameter  interval  to  three-dimensional  space 

r(u)  =  [x(u)   y(u)   z{u)\T . 

The  points  on  a  curve  are  regular  as  long  as  at  least  one  of  the  first  derivatives  is  not  equal 
to  zero.  In  other  words,  a  point  is  singular  if  dx/du  =  dy/du  =  dz/du  =  0.  The  curve  can 
be  reparameterized  if  u  can  be  expressed  as  u  =  f(u\)  as  long  as  du/du\  ^  0. 

Arc  length  is  defined  as  the  distance  along  a  curve  between  parameter  values  uo  and  u, 
and  can  be  represented  as 

s{u)  =        y/x2  +  y2  +  z2  du=   /    Vr~~tdu.  (2.1) 


In  this  thesis,  derivatives  with  respect  to  the  arc  length  s  will  be  represented  by  a  prime 
and  derivatives  with  respect  to  u  will  be  represented  by  a  dot.  Derivatives  of  arc  length  s 
with  respect  to  parameter  u  and  vice  versa  are 


ds  r^^ 

s    =     — -  =  vr  •  r  =  r 


du 

ds  TV 


s     = 


u' 


du       y/r  ■  r 
du        1 
~ds~  ==  \F\ 
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u"     = 


duf_ 

ds 


r  •  r 

(r-r): 


The  unit  tangent  vector  t  to  a  curve  is  defined  by  the  unit  vector  that  passes. through 
two  points  in  the  curve,  r(s)  and  r(s  +  h)  (see  Figure  2-1),  as  h  approaches  zero;  i.e. 

r(s  +  h)  -r(s)       dr        , 

t     =     lim =  —  =  r  (5) 

ft-*0  h  ds 

dr(u)  du        t{u) 


du    ds       |r(u)| 

This  final  representation  clearly  shows  that  t  is  a  unit  vector.  The  line  passing  through  the 
point  r(un)  in  the  direction  of  t(un)  is  the  tangent  to  the  curve  at  un. 


tangent  line. 


Figure  2-1:  Tangent  Vector 


All  vectors  normal  to  t(un)  at  un  lie  in  a  plane  called  the  normal  plane.  Much  as 
the  tangent  is  defined  by  two  points  that  approach  one  another  on  a  curve,  the  osculating 
plane  is  defined  by  three  points  that  approach  one  another.  In  other  words,  the  osculating 
plane  is  the  plane  in  which  the  curve  lies  at  a  point  on  the  curve.  If  the  curve  is  planar, 
the  entire  curve  lies  in  a  single  osculating  plane,  and  the  osculating  plane  is  constant  for 
the  entire  curve.  To  mathematically  define  the  osculating  plane  following  the  derivation  in 
Kreyszig  [26,  p. 31],  take  three  points  at  parameter  values  u,  u  +  h\  and  u  +  h.2-  The  plane 
can  be  defined  by  the  two  linearly  independent  vectors  ai  =  r(u  +  hi)  —  r(u),  i  =  1,2  or 
linear  combinations  thereof.  Let 


VW  = 


v(u  +  hi)  -  t(u) 


and 


w  = 


2(v 


(2)  _  v(l) 


/l2  —  ^1 
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The  Taylor  series  expansion  of  r(u  +  hi)  is 


r(u  +  hi)  =  r(u)  +  /izr(u)  +  -^-r(t*)  +  o(h-) 
where  o(hn)  is  a  vector  of  Landau  symbols  o(hn)  with  the  property  that  [41,  p.3] 

Urn  °m  =  0. 


/i->o    hn 


(2.2) 


Using  (2.2),  we  obtain 


r(i)     = 


r(u  +  h\)  —  r(u) 


r(u)  +  hir(u)  +  j/iff (u)  +  o(h\)  -  r(u) 


and 


w 


=    r(u)  + -hir(u)  + -^ 


2(V(2)_V(D) 


/l2  —  /ll 

2 


ho  —  h 


o{hl) 


r(u)  +  -/i2f  (u)  + 

2  Al2 


=     f(u) 


2(/iio(/ij)-M(^)) 
^1^2(^2  —  /il) 


*(«)  +  \hxr(u)  +  2iMl 


As  /ij  — >  0,  v^1)  ->  t(u)  and  w  — >  r(it).  Therefore,  the  osculating  plane  is  defined  by  r(u) 
and  r(u).  Note  that  r(u)  must  be  at  least  twice  continuously  differentiable. 


osculating  plane 


Figure  2-2:  Osculating  Plane 

The  principal  normal  is  defined  as  the  intersection  of  the  osculating  and  normal  planes. 
Using  an  arc  length  parameterization,  we  find  that  differentiating  t  •  t  =  1  with  respect  to 
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s  yields 

t  ■  t'  =  0, 

so  t'  is  perpendicular  to  t  and  lies  in  the  normal  plane.  Since  t  =  r'  =  rit', 

t'  =  ru"  +  r(u')2. 

Thus  t'  lies  in  the  plane  of  r  and  r,  and  therefore  in  the  osculating  plane.  Since  t'  lies  in 
both  the  osculating  and  normal  planes,  a  unit  vector  in  the  direction  of  t'  is  the  unit  normal 
vector 


t' 


n  = 


The  magnitude  of  t'  is  called  the  curvature  and  is  represented  as  k.  The  curvature  vector 
k  is  represented  as 


k  =  t'  =  nn. 


(2.3) 


A  physical  sense  of  the  curvature  can  be  found  from  the  simple  example  of  a  line,  which 
has  a  constant  tangent  vector,  so  t'  is  zero.  Therefore,  k  =  |t'|  =  0,  and  a  line  has  zero 
curvature.  Parameterized  by  u,  k  and  «  can  be  represented  as  [41,  p.  16] 


*■>      ;.  n  ,«,  '\2       (r  •  r)f  -  (f  •  f)r 

t  —  ru    +  r[u  )    =  - 


—     («n)  •  (/in)  = 


(r-r)2 
(r  •  r)r  —  (r  •  r)r 


'(r  •  r)r  —  (r  •  r)r 


(r-r): 


(r  x  f )  ■  (r  x  r) 
(r-r)3 


(r-r)2 

where  the  Lagrange  identity  (a  x  b)  •  (a  x  b)  —  (a  •  a)(b  •  b)  -  (a  ■  b)2  is  used. 

The  normal  to  the  osculating  plane  is  the  binormal,  which  is  defined  as 

b  =  t  x  n. 

Since  these  three  vectors  are  unit  vectors  that  are  mutually  orthogonal,  they  satisfy  the 
relations 


b   b  =  1 

n    n  =  1 

t-t  =  1 

bn  =  0 

nt  =  0 

t  •  b  =  0. 

Differentiating  b  •  b  =  1  and  b  •  t  =  0  yields 

b'  •  b  =  0 
and 

b'  •  t  =  -b  •  t'  =  -b  •  «n  =  -«(b   n)  =  0. 
Therefore,  b'  is  orthogonal  to  both  b  and  t  and  must  then  be  parallel  to  n.  Let 

b'  =  -rn 
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where  r  is  the  torsion  of  the  curve.  A  physical  sense  of  torsion  can  be  found  from  the  fact 
that  a  curve  with  zero  torsion  and  nonzero  curvature  is  a  planar  curve. 

Finally,  an  equation  for  the  derivative  of  the  normal  vector  corresponding  to  the  rate  of 
change  of  the  tangent  plane  [26,  p. 41]  can  be  determined.  Differentiating  n  •  n  =  1  yields 
n  •  n'  =  0.  Therefore,  n'  is  orthogonal  to  n  and  must  satisfy  the  relationship 

n'  =  at  +  /3b.  (2.4) 

Multiplying  equation  (2.4)  by  t  and  b  respectively  yields 

a  =  n'  •  t     and     j3  =  n'  •  b. 

Differentiating  n  •  t  =  0  yields 

a  =  n'  •  t  =  — n  •  t'  =  — n  •  nn  =  —k. 

Similarly,  differentiating  n  •  b  =  0  yields 

(3  =  n'  •  b  =  — n  •  b'  =  — n  •  (— rn)  =  r. 

Thus, 

n'  =  —  nt  +  rb. 

The  representations  of  the  derivatives  of  the  three  defining  vectors  for  a  curve  are  termed 
the  Serret-Frenet  formulae  [41,  p. 18] 

t'       =       Kli 

b'     =     -rn 

n'     =     -Kt  +  rb, 

which  describe  all  aspects  of  a  space  curve,  as  shown  in  Figure  2-3. 


Osculating  Plane 


Normal  Plane 


Tangent  Plane 


Figure  2-3:  Space  Curve 
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2.3     Basic  Theory  of  Surfaces 


A  general  parametric  surface  can  be  defined  as  a  vector-valued  mapping  from  a  two- 
dimensional  parametric  uv-spa.ce  to  three-dimensional  space 

R(u,v)  =  [x(u,v)  y(u,v)  z(u,v)}   , 

where  the  two  parameters  u  and  v  lie  in  the  closed  intervals  u\  <  u  <  U2  and  v\  <  v  <  V2- 
The  surface  is  regular  at  all  points  where  the  matrix 


M  = 


dx  djl  dz_ 

du  du  du 

dx  chi  ch 

dv  dv  dv 


(2.5) 


has  rank  2,  or  equivalently  R^  x  R„  ^  0  [41,  p. 56].  Partial  derivatives  of  a  surface  will  be 
represented  as 

dR  OR 

Ru  =  —    and    R^  =  — . 
au  ov 

A  curve  on  the  surface  is  represented  as  r(u,  v)  with  u  =  u(t)  and  v  =  v(t)  or  equivalently 
R(u(t),  v(t))  =  r(t).  This  occurs  when  the  rank  of  the  matrix  (2.5)  is  1  at  every  point. 
When  we  keep  u  constant  in  R(u,v)  by  setting  u  =  un,  we  have  an  equation  for  a  curve 
that  depends  only  on  v,  and  is  thus  called  an  isoparameter  curve.  Similarly,  v  =  vn  is  also 
an  isoparameter  curve. 

The  vector 

is  tangent  to  the  curve  when  u  and  v  are  functions  of  t.  This  is  of  course  also  a  tangent  to 
the  surface.  The  tangent  plane  to  the  surface  at  any  point  is  determined  by  the  two  vectors 
Ru  and  R^.  The  normal  to  the  surface  is  orthogonal  to  this  plane.  A  unit  vector  in  the 
direction  of  the  surface  normal  is  given  by 

N  ■  iirrsn-  (2J) 

|xvu  x  rt^l 

Note  that  this  definition  of  the  normal  vector  is  undefined  when  R^  x  Ry  =  0,  supporting 
the  assertion  that  a  point  on  the  surface  is  singular  if  R„  x  R^  =  0.  This  may  be  due  to 
the  shape  of  the  surface  or  the  choice  of  parameterization.  In  Section  5.3.5,  an  alternative 
definition  of  N  is  explored  for  degenerate  points. 

Using  the  definition  of  arc  length  (2.1),  the  distance  between  two  points  on  a  curve  on 
the  surface  is  found  by  integrating 


ds  =  y/dxdx  +  dydy  +  dzdz  =  VdR  ■  dR.  (2.8) 

Combining  (2.6)  and  (2.8)  yields  the  first  fundamental  form  of  the  surface 

I  =  ds2  =  {Rudu  +  Rydv)  ■  (Rudu  +  Rydv)  =  Edu2  +  2Fdudv  +  Gdv2 
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where 

E  =  Ry  •  Ru,      F  =  Ru  •  Rt,,      G  =  R^  ■  R^. 

Note  that 

EG  —  F      =     (Ru  •  RU)(R1,  ■  Rt;)  —  (R^  •  Ry) 


=     (Ku  x  Rv)  •  (Ru  x  Ry) 

=     KRuXR,)!2.  (2.9) 


Thus,  EG  —  F2  is  always  positive  if  the  surface  is  regular. 


The  first  fundamental  form  enables  us  to  determine  arc  lengths,  angles  and  areas  on  the 
surface.  In  order  to  determine  curvatures,  the  second  fundamental  form  is  required.  First, 
the  curvature  vector  of  the  curve  on  the  surface,  found  from  equation  (2.3),  is  split  into  its 
normal  and  tangential  components 


ds 


K  —     t     —  iCfj  -+■  Kg  —  ACjjiN  t  K^ 


where  kn  is  the  normal  curvature  vector  and  Kn  is  the  normal  curvature  of  the  curve  on 
the  surface  at  that  point.  The  tangential  portion  of  the  curvature  vector,  k5,  is  called  the 
geodesic  curvature  vector.  It  will  be  discussed  in  more  detail  later. 


If  the  angle  between  the  surface  normal  and  the  normal  to  the  curve  is  designated  7  as 
shown  in  Figure  2-4,  then  [26,  p.  118] 

n  •  N  =  cos(7)|n||N|  =  cos(7) 

from  the  definition  of  a  dot  product.  Since  t'  =  R"  =  k  ■  n, 

KCOs(-y)     =  k  ■  n  •  N 

=  R"-N 

=  (Rudu  +  Rydv)'  ■  N 

=  (Ruududu  +  2Ruvdudv +  Rvvdvdv  +  Rud2u  +  Rvd2v) -N.       (2.10) 

Since  Ry  •  N  =  R^  •  N  =  0,  equation  (2.10)  reduces  to  the  second  fundamental  form, 
represented  by 

77    =     (Ruududu  +  2Ruvdudv  +  Ryydvdv)  ■  N 
=    Ldudu  +  2Mdudv  +  Ndvdv 


where 


L    =    Ruu-N  (2.11) 

M    =    RuvN       and  (2.12) 

N    =    R^-N.  (2.13) 
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Also  note  that  since  R^  •  N  =  0, 

<J(R^  •  N) 
du 


=    KuU  ■  N  4-  Ru  •  Nu  =  0, 


and  therefore 

Similarly,  R^  •  N  =  -Ry  ■  Nv  =  -R^  •  Nu  and  R^  ■  N  =  -R,;  •  N„.  Therefore, 

L    =     -Ku-Nu 
M    =    -Ry  •  N,,  =  -R^,  •  Nu 

=    --(Ru-Nv  +  R^-Nu)       and 

N    =    -R„-N„. 


(2.14) 
(2.15) 

(2.16) 


N 


Figure  2-4:  Curvature  Vectors 

To  find  the  normal  curvature,  we  begin  with  the  fact  that  the  surface  normal  is  orthog- 
onal to  the  tangent  to  the  curve.  Differentiating  t  •  N  =  0  with  respect  to  s  yields 

dt  dN  _ 

ds  ds 


Thus  we  have 


k    IV-  —    1ST-  -      —  dR     dN 

ds  ds  ds      ds 


Therefore, 


_  dR    dN       (Kudu  +  Rydv)  ■  (Nudu  +  Nvdv)  II 

ds      ds  ds2  I 

Ldudu  +  2Mdudv  +  Ndvdv 

Edudu  +  2Fdudv  +  Gdvdv 
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L  +  2MX  +  NX2 

E  +  2FX  +  GX2  [  '     ' 

where  A  =  dv/du.   The  negative  sign  in  equation  (2.17)  ensures  that  if  Kn  is  positive,  the 
center  of  curvature  lies  opposite  to  the  direction  of  the  surface  normal. 

Any  point  on  the  surface  has  many  curves  passing  through  it,  each  of  which  has  an 
associated  Kn.  The  maximum  and  minimum  values  of  Kn  at  a  point  are  called  the  principal 
curvatures  of  the  surface  at  that  point  and  are  designated  K\  and  «2-  The  maximum  and 
minimum  values  of  nn  occur  when  ^  =0,  [11,  p.  112] 

dK 

-f  =  (£  +  2FA  +  GX2){2M  +  2NX)  -  (L  +  2MX  +  NX2)(2F  +  2GX)  =  0.        (2.18) 
dX 

Using  (2.18)  in  (2.17)  yields 

L  +  MX         M  +  NX 


Kn  E  +  FX~'=      F  +  GX 


given 


{E  +  2FX  +  GX2)     =     {E  +  FX)  +  X(F  +  GX)    and 
(L  +  2MA  +  JVA2)     =     (L  +  MX)  +  X(M  +  NX). 

Therefore  the  extreme  values  of  Kn  satisfy  the  two  simultaneous  equations 

(L  +  KnE)du  +  (M  +  KnF)dv    =    0      and  (2.19) 

(M  +  KnF)du  +  {N  +  KnG)dv    =    0.  (2.20) 


These  equations  can  be  simultaneously  satisfied  if  and  only  if 

=  0. 


L  +  KnE     M  +  KnF 
M  +  KnF     N  +  KnG 


The  discriminant  of  this  quadratic  equation  in  Kn  is  greater  than  or  equal  to  zero.  Therefore 
the  equation  has  either  two  distinct  roots  Kmax  and  Kmin,  the  maximum  and  minimum 
principal  curvatures,  or  double  roots  «n,  the  normal  curvature  at  an  umbilical  point.  The 
corresponding  directions  A  define  directions  in  the  wu-plane.  The  corresponding  directions 
in  the  tangent  plane  are  called  principal  directions  of  curvature  and  are  always  orthogonal 
except  at  umbilical  points,  where  nmax  and  Kmm  are  identical.  In  the  special  case  where  the 
identical  principal  curvatures  vanish,  the  surface  becomes  locally  flat.  The  two  roots  are 
given  by 


>^max    =    H  +  \/H2  -K      and  (2.21) 

Kmin     =    H  -  \/H2-K  (2.22) 

where  K  is  the  Gaussian  curvature  and  H  is  the  mean  curvature  defined  by 

LN  -  M2 

K  =  ^rrp  <2-23> 

2FM  -EN-GL  ,„„„. 

H    m         2{EG-F>)      ■  (224) 
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Using  (2.21)  and  (2.22)  in  (2.23)  and  (2.24),  it  can  easily  be  shown  that 


maxnmin 


K      =      K 

H      =      -(Kmax+Krmn). 

Note  that  at  an  umbilical  point,  H2  =  K  and  at  a  flat  point,  K  =  H  =  0. 

The  normal  curvature  can  be  represented  in  terms  of  the  principal  curvatures  «i  and  «2 
as 

Kn  =  K\  cos2 (a)  +  K2  sin2 (a)  (2.25) 

where  a  is  the  angle  between  the  direction  dv/du  and  the  direction  dv  =  0.  This  is  known 
as  Euler's  theorem  [26]. 

Recall  that  the  curvature  vector  of  a  curve  on  the  surface  can  be  split  into  its  normal 
and  tangential  components 

K  =  Kfj  T"  Kg, 

where  the  tangential  component,  ks,  is  known  as  the  geodesic  curvature  vector.  Let  u  be 
a  unit  vector  in  the  tangent  plane  perpendicular  to  t  such  that  t,  u  and  N  form  a  right- 
handed  coordinate  system  in  that  order  as  shown  in  Figure  2-4.  The  geodesic  curvature  ng 
can  then  be  found  from 

Kg    =    KgU. 

Additionally,  since  the  magnitude  of  «n  can  be  represented  as  |kcos7|,  the  magnitude  of 
Kg  can  be  represented  as  |«sin7|. 

Geodesic  curvature  can  be  represented  using  only  the  first  fundamental  form.    Since 
u  •  N  =  0  and  u  •  u  =  1, 


dt 
ds 


Kj  =  u-k  =  u--  =  u-t'  =  (Nxt)-t'=  (tt'N).  (2.26) 


The  unit  tangent  vector  of  the  curve  on  the  surface  is  given  by 

t  =  TLuU  +  R^t/; 

thus 

t'  =  Huu(u')2  +  2RuvtiV  +  RvV(v')2  +  Kuu"  +  Kyv". 

Therefore  equation  (2.26)  can  be  rewritten  as 

Kg     =     (t  x  t')  •  N 

=     [(Rut*'  +  R^v')  x  {Kuuiu')2  +  2RuvuV  +  IW*/)2  +  R^"  +  Rvv")]  ■  N 
=     [(R,,  x  Ruu)(u')3  +  (2Ru  xKuy+KyX  RuU)(u')2v' 

+  (Ru  xR,„+2R„x  KuV)u'(v')2 

+  (Rr  x  RyvUv')3  +  (Kux  R„)(uV'  -  uV)]  ■  N.  (2.27) 
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Each  of  the  coefficients  of  the  combinations  of  u  and  v  derivatives  in  (2.27)  can  be  written 
in  terms  of  the  coefficients  of  the  first  fundamental  form  and  their  derivatives.  For  example, 

(R.xR_).N    -  ^eg  _  ^ 

(Ru  •  R^j^R^u  •  Rt,)  —  (Ftu  *  R-uu)(Ru  "  *M;) 


■R.) 

VEG 
-(Ru 

-F2 

•  Ruu)-? 

2(£(RUW 

7)V^G- 

,F2 

{2EFU  - 

EEU 

2(£G- 

-F2) 

v/FG- 

,F2 

2(EG  -  F2) 


=  t2uVeg-f2 

where  r2:  is  one  of  the  Christoffel  symbols  defined  by  the  coefficients  of  the  first  fundamental 
form  E,  F  and  G  and  their  derivatives,  Eu,  Fu,  Gu,  Ev,  Fv  and  Gv 

i  GEU  —  2FFU  +  FEV  2        2EFU  —  EEV  +  FEU 

11  ~           2(EG  -  F2)        '  "  ~         2{EG-F2) 
ri             GEV  —  FGU  2        EGU  —  FEV 

12  "     2{EG-F2Y  12  ~  2{EG-F2) 


,!  2GFV  —  GGU  +  FGV  2        EGV  —  2FFV  +  FGU 

22  2{EG  -  F2)        '  22  ~         2(EG  -  F2) 


Similarly  reducing  all  the  other  coefficients  in  (2.27)  yields 

«*   =   [r21(nr  +  (2r22-r!j(ur^'  +  (r22-2r;2KK)2-rLK)3 

+u'v"  -  u"v']  s/EG  -  F2.  (2.28) 

A  geodesic  is  defined  as  a  curve  with  zero  geodesic  curvature  [41].  Straight  lines  on  a 
surface  are  geodesies  since  the  curvature  vector  k  vanishes.  For  geodesies  that  are  curved, 
the  curvature  vector  coincides  with  the  surface  normal  vector  and,  since  the  curvature  vector 
lies  in  the  osculating  plane,  the  osculating  and  tangent  planes  are  normal.  The  equation  for 
the  geodesic  can  be  obtained  by  setting  ng  equal  to  zero  in  equation  (2.28),  which  yields 

u'V  -  u'v"  =  r2n(u')*  +  (2r22  -  rh)(u')V  +  (r22  -  2r}»')2  -  r\2(v'f 

assuming  that  the  geodesic  is  everywhere  regular  so  EG  —  F2  is  always  positive.  An  al- 
ternative representation  for  a  geodesic  is  given  by  the  set  of  coupled  second  order  ordinary 
differential  equations  [41] 

d2u      „■,    ( du\2      „„n   dudv      „■,    f  dv\2  . 

^  +  r"U)  +2r«55  +  r»U)  =°  <2'29) 

^  +  r?1(^2  +  2r?2^  +  r|2f^2  =  o  (2.30) 


ds2         n  \dsj  12dsds        ll  \ds 

where  the  two  equations  are  related  by  the  condition  ds2  =  Edu2  +  2Fdudv  +  Gdv 
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Chapter  3 

Curve  and  Surface  Representation 


3.1     Introduction 

This  chapter  reviews  the  representation  of  curves  and  surfaces  in  Bezier  and  B-spline  forms 
and  treats  the  special  properties  associated  with  each.  In  addition,  algorithms  for  knot 
insertion  and  removal,  curve  evaluation  and  splitting  are  summarized.  The  descriptions  are 
based  on  texts  by  Hoschek  and  Lasser  [22],  Piegl  and  Tiller  [36]  and  Yamaguchi  [44]. 


3.2     Bezier  Curves  and  Surfaces 

Bernstein  Polynomials 

The  Bernstein  polynomials  are  defined  as 

5,,n(u)  =        "^,(1  -  u)n~lu\        i  =  0, . . . ,  n,  (3.1) 

and  have  several  properties  of  interest.  The  property  of  positivity  states  that  each  polynomial 
factor  is  non-negative  such  that 

Bi,n(u)  >  0,  0  <  u  <  1 

for  all  i  and  n.  The  partition  of  unity  property  states  that  the  Bernstein  polynomials  sum 
to  1  for  all  0  <  u  <  1,  or 

n 

£^,n(u)  =  l. 
i=0 

The  derivative  of  a  Bernstein  polynomial  is 

dBj,n(u) 

— — —  =  n^-i^-itu)  -  Bltn-i{u)\ 

where  B_iiTl_i  =  B„)n_i  =  0.  The  linear  precision  property  of  a  Bernstein  polynomial 


n 
I 


U  =  V  -Bin(u) 


t=on 
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allows  the  monomial  u  to  be  expressed  as  the  weighted  sum  of  Bernstein  polynomials  with 
coefficients  evenly  spaced  in  the  interval  [0,1]. 

Bezier  Curves 

A  Bezier  curve  is  a  spline  curve  that  uses  the  Bernstein  polynomials  as  a  basis.  A  Bezier 
curve  of  degree  n  (order  n  +  1)  is  represented  by 


r(u)  =  £btBt,n(ti),        0<«<1. 


t=0 

The  coefficients,  b2,  are  the  control  points  that  determine  the  shape  of  the  curve.  Lines 
drawn  between  consecutive  control  points  of  the  curve  form  the  control  polygon.  A  cubic 
Bezier  curve  is  shown  in  Figure  3-1.  Bezier  curves  have  the  following  properties: 

•  The  first  and  last  control  points  are  the  endpoints  of  the  curve.    In  other  words, 
b0  =  r(0)  and  bn  =  r(l). 

•  The  curve  is  tangent  to  the  control  polygon  at  the  endpoints. 

•  The  convex  hull  property  states  that  the  entire  curve  is  contained  within  the  convex 
hull  of  the  control  points. 

•  The  first  derivative  of  a  Bezier  curve  is  represented  by 
r(u)  =  ^^  =  n  ]T  (bi+i  -  bi)Bijn-i{u),       0  <  u  <  1. 


du 


i=0 


bi. 


Convex  Hull 


Control  Polygon 


b2 


Figure  3-1:  A  Cubic  Bezier  Curve  with  Control  Polygon 

A  Bezier  curve  can  be  evaluated  at  a  specific  parameter  value  uq  and  the  curve  can  be 
split  at  that  value  using  the  deCasteljau  algorithm,  where  [22,  p. 125] 


.*-i 


bf(w0)  =  (1  -tio)b$_{  +u0b*   \     fc  =  l,2,...,n,     i  =  k,...,n 
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is  applied  recursively  to  obtain  the  new  control  points.    The  algorithm  is  illustrated  in 
Figure  3-2,  and  has  the  following  properties: 

•  The  values  b-*  are  the  original  control  points  of  the  curve. 

•  The  value  of  the  curve  at  parameter  value  ito  is  b™. 

•  The  curve  can  be  represented  as  two  curves,  with  control  points  (bj ,  b", . . . ,  b£)  and 
(bn    bn_1  b°) 

•  Since  this  is  merely  a  change  in  the  basis  representation,  the  curve  itself  remains 
parametrically  and  geometrically  unchanged. 


b?=bf=b]; 


b§=b§=bj=bg 


Figure  3-2:  The  deCasteljau  Algorithm 


Bezier  Surfaces 


A  tensor  product  surface  is  formed  by  moving  a  curve  through  space  while  allowing  defor- 
mations in  that  curve.  This  can  be  thought  of  as  allowing  each  control  point  b;  to  sweep  a 
curve  in  space.  If  this  surface  is  represented  using  Bernstein  polynomials,  a  Bezier  surface 
is  formed,  with  the  following  formula: 

m      n 

R(u,  v)  =  J3 5Z  by'Bi,m(«)5j,f.(«),  0<U,V<1. 

i-Oj-0 

Here,  the  set  of  lines  drawn  between  consecutive  control  points  bij  is  referred  to  as  the 
control  net.  An  example  of  a  bi-quadratic  Bezier  surface  with  its  control  net  can  be  seen  in 
Figure  3-3.  The  following  conditions  apply: 

•  The  boundary  isoparametric  curves  (u  =  0,  u  =  1,  v  =  0  and  v  =  1)  have  the  same 
control  points  as  the  corresponding  boundary  points  on  the  net. 

•  The  corners  of  the  surface  coincide  with  the  corner  points  of  the  net,  and  the  deriva- 
tives in  the  u  and  v  directions  at  the  corners  are  tangent  to  the  net. 

•  Bezier  surfaces  exhibit  the  convex  hull  property. 
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b'oi 
Figure  3-3:  A  Bezier  Surface  with  Control  Net 

Continuity  Conditions 

Bezier  curves  can  represent  complex  curves  by  increasing  the  degree  and  thus  the  number 
of  control  points.  Alternatively,  complex  curves  can  be  represented  using  composite  curves, 
which  can  be  formed  by  joining  several  Bezier  curves  end  to  end.  If  this  method  is  adopted, 
the  continuity  between  consecutive  curves  must  be  addressed. 

One  set  of  continuity  conditions  are  the  geometric  continuity  conditions,  designated  by 
the  letter  G  with  an  integer  exponent.  Position  continuity,  or  G°  continuity,  requires  the 
endpoints  of  the  two  curves  to  coincide, 

rW(l)  =  r<2>(0). 

The  superscripts  denote  the  first  and  second  curves.  Tangent  continuity,  or  G1  continuity, 
requires  the  tangents  of  the  curves  to  be  in  the  same  direction, 

r(1)(l)=ait 
r(2)(0)  =a2t 

where  t  is  the  common  unit  tangent  vector.  Gl  continuity  is  important  in  minimizing  stress 
concentrations  and  preventing  flow  separation. 

Curvature  continuity,  or  G2  continuity,  requires  the  center  of  curvature  to  move  contin- 
uously past  the  connection  point, 

r(2)(0)  =  ^)2r(^)(l)+^1)(l). 

G2  continuity  is  important  for  aesthetic  reasons  and  for  preventing  fluid  flow  separation. 

More  stringent  continuity  conditions  are  the  parametric  continuity  conditions,  where 
Ck  continuity  requires  the  A;th  derivative  of  each  curve  to  be  equal  at  the  joining  point.  In 
other  words, 

dkrM(l)  __  rfM2>(0) 
duk  duk 
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The  Cl  and  C2  continuity  conditions  for  consecutive  segments  of  a  composite  degree  n 
Bezier  curve  can  be  stated  as  [44,  p.  215] 


hi+i  (bn»  ~  bnj_i)     =    ht  {bni+i  -  bm) ,       and 


(3.2) 


hi+] 


bni-1  +  -r — (bni-i  -  bni-2)      =      bnj+i  + (b„j+i  -  bnj+2)  (3.3) 

where,  for  the  ith  Bezier  curve  segment  with  u  values  between  knot  values  tx  and  ti+\ 
inclusively,  ht  —  t{+i  —  t%.  See  Figure  3-4. 


bni-2- 


•■--..hi 


'ni-3 


Figure  3-4:  Continuity  Conditions 


3.3     B-Spline  Curves  and  Surfaces 

The  Bezier  representation  has  two  main  disadvantages.  First,  the  number  of  control  points 
is  directly  related  to  the  degree.  Therefore,  to  increase  the  complexity  of  the  shape  of  the 
curve  by  adding  control  points  requires  increasing  the  degree  of  the  curve  or  satisfying  the 
continuity  conditions  between  consecutive  segments  of  a  composite  curve.  Second,  changing 
any  one  control  point  affects  the  entire  curve  or  surface,  making  design  of  specific  sections 
very  difficult.  These  disadvantages  are  remedied  with  the  introduction  of  the  B-spline 
representation. 

A  B-spline  (basis-spline)  curve  is  a  piecewise  polynomial  curve  meaning  it  is  made  up  of 
polynomial  curve  segments.  It  is  a  spline  curve  with  a  different  basis  than  the  Bezier  curves, 
although  a  Bezier  curve  is  a  special  case  of  a  B-spline  curve.  The  B-spline  representation 
has  the  advantage  that  the  control  points  affect  the  curve  only  locally,  and  a  B-spline  of 
order  k  may  have  as  many  control  points  as  required  to  describe  the  curve. 

A  B-spline  curve  is  a  series  of  polynomial  segments  of  degree  k  —  1  joined  together  at 
knots  tk.  A  B-spline  of  order  k  (degree  k  —  1)  with  n  +  1  control  points  has  the  equation 

n 

r(u)  =  ^Pi^W^),        n>k-l,        ue[<w»Wl] 

i-0 
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where  tk  is  a  member  of  the  knot  vector  with  n  +  k  +  1  elements 

T  =  (^0^1)  •  •  ■  i£fc-l)£jfci*Jfc+li  •  •  •  ,tn-i,tn,tn+\,  •  •  •  ,tn+h),  to  <  t\  <  .  .  .  <  tn+k- 

The  basis  functions,  Nijk(u),  are  given  by: 

1     for  tx  <  u  <  U+i 


for  k  =  1,  and 


Ni,i(u)  -  . 


NijkM  =      "     V  Ni,k-i(u)  +    *i+k   4U  Ni+l>k^(u) 

H+k-l  ~  H  H+k.  ~  H+l 


for  k  >  1  and  i  =  0, 1, ...  ,n.  These  equations  have  the  following  properties  [22,  p.  168]: 

•  Nitk[u)  >  0,  for  U  <  u  <  tl+k 

•  Nitk{u)  =  0,  for  to  <u<tt,  and  tl+k  <u<  tn+k 

•  £"=o-Nifc(u)  =  1,  for  u  e  [tk-i,tn+1],  and 

•  Nlk(u)  has  continuity  Ck~2  at  each  simple  knot  £/. 

The  first  three  properties  show  positivity  and  partition  of  unity  for  the  B-Spline  basis 
functions. 


Pr-3 


Pr-1 


Figure  3-5:  A  Cubic  B-spline  Curve  Segment  with  Control  Polygon 

Local  control  means  that  a  single  segment  of  the  B-spline  curve  is  controlled  only  by  the 
nearest  n  points,  and  that  any  control  point  affects  the  nearest  n  spans.  In  other  words, 
changing  pt  affects  the  curve  in  the  parameter  range  U  <  u  <  tl+k  and  the  curve  at  a  point 
u  where  tr  <  u  <  tr+\  is  determined  completely  by  the  control  points  pr_^_!), . . . ,  pr. 

The  convex  hull  property  for  B-splines  applies  locally,  so  that  a  span  lies  within  the 
convex  hull  of  the  control  points  that  affect  it.  This  provides  a  tighter  convex  hull  property 
than  that  of  a  Bezier  curve,  as  can  be  seen  in  Figure  3-6. 
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Figure  3-6:  Convex  Hull  Property  for  a  Cubic  B-Spline  Curve 


Increasing  the  multiplicity  of  a  knot  reduces  the  continuity  of  the  curve  at  that  knot. 
Specifically,  the  curve  is  (k—p—1)  times  continuously  differentiable  at  a  knot  with  multiplic- 
ity p  (<  k),  and  thus  has  C^k~p~1^  continuity.  Therefore,  the  control  polygon  will  coincide 
with  the  curve  at  a  knot  of  multiplicity  k  —  1,  and  a  knot  with  multiplicity  k  indicates  C_1 
continuity,  or  a  discontinuous  curve. 

As  you  can  see  in  Figure  3-5,  the  endpoints  of  the  curve  may  not  coincide  with  the 
control  polygon.  However,  repeating  the  knots  at  the  end  k  times  will  force  the  endpoints 
to  coincide  with  the  control  polygon.  Thus  a  knot  vector  described  by 

T  =  (£o,£i, . . .  ,iA.—i,  tk,tk+i, . . .  ,tn-\,tn,  tn+\, . . . , £n+jfc) 

V v '     * v '     V v ' 

k  equal  knots    n  —  k  +  1  internal  knots  k  equal  knots 

will  have  to  —  t\  =  . . .  =  t^-i  and  tn+\  =  tn+2  =  •  •  •  =  tn+k  and  hence  the  first  and  last 
control  points  of  the  curve  coincide  with  the  endpoints  of  the  curve. 

From  this  discussion,  it  can  be  seen  that  a  Bezier  curve  of  order  k  (degree  k  —  1)  is  a 
B-spline  curve  with  no  internal  knots  and  the  end  knots  repeated  k  times.  The  knot  vector 
is  thus 

T  =  (^Oi^l,  •  •  •  ,£fc-l,*n+l>  •  •  •  ,tn+k) 
«. „ '  > w ' 

k  equal  knots      k  equal  knots 
where  n  +  k  +  l=  2k  or  n  =  k  —  1. 

A  B-spline  surface  is  a  tensor  product  surface  using  a  B-spline  basis.  This  is  represented 
mathematically  as 

m     n 

R(w,  v)  =  J2  J]  PijNiyk(u)Nj}i(v). 

i=0j=0 
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Knot  Insertion 

A  knot  can  be  inserted  into  a  B-spline  curve  without  changing  the  geometry  of  the  curve. 
The  new  curve  is  identical  to  the  old  one,  with  a  new  basis  where 


n+l 


^2piNi)k(u)     becomes     ^p^^w) 

i-0  i-0 

over  T  =  [t0,  *i,  —  ,*i,  *i+r,  •  •  •]  over  T  =  [t0,  tu 


,ti,t,ti+i, . . .] 


when  a  new  knot  i  is  inserted  between  knots  t\  and  ti+\.  The  deBoor  algorithm  is  used  to 
insert  the  knot,  and  is  represented  as 


a = a  -  c^pt; + cid-1    i>i-k+2 


(3.4) 


where 


Qi 


"o-ti 


tl+k-l—U 


i<l-k+l 
l  +  l  <i 

l-k+2<i<l 


To  evaluate  a  B-spline  curve  at  a  specific  parametric  value  or  to  split  the  curve  at  that 
value,  the  deBoor  algorithm  is  used  to  repeatedly  insert  the  same  knot  until  the  control 
polygon  coincides  with  the  curve  at  the  parametric  value  wo,  or  in  other  words,  until  the 
multiplicity  of  the  knot  at  «o  =  k  —  1.  Then, 

P?     =    Pi 

PfclJ     =    r(u0). 

A  B-spline  curve  is  C°°  continuous  in  the  interior  of  a  span.  Inserting  a  knot  does  not 
change  the  curve,  so  it  does  not  change  the  continuity.  However,  if  any  of  the  control  points 
are  moved  after  knot  insertion,  the  continuity  at  the  knot  will  become  Ck~p~l,  where  p  is 
the  multiplicity  of  the  knot.  Figure  3-7  illustrates  a  single  insertion  of  a  knot  at  parameter 
value  wq,  resulting  in  a  knot  with  multiplicity  of  one. 


P? 


P^ 


,p£ 


Po=P8 


P°3=Pl 


Figure  3-7:  The  deBoor  Algorithm 
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The  B-spline  curve  can  be  subdivided  into  Bezier  segments  by  knot  insertion  at  each 
internal  knot  until  the  multiplicity  of  each  internal  knot  is  the  degree  k.  The  above  algorithm 
is  also  known  as  Boehm's  algorithm  [4,  5]. 

Knot  Removal 

Knot  removal  is  the  reverse  process  of  knot  insertion.  This  thesis  uses  the  knot  removal 
algorithm  developed  by  Tiller  [42],  in  which  a  knot  is  removed  if  and  only  if  it  leaves  the 
curve  or  surface  geometrically  and  parametrically  unchanged. 

To  demonstrate  the  process,  this  example  uses  a  cubic  B-spline  curve  r(u)  given  by 
control  points  (pg,  ■ . . ,  p°)  and  knot  vector  (£o,  •  •  •  *io)  where  to  =  . . .  =  t3  =  0,  t±  —  t$  = 
<6  =  1  and  £7  =  . . .  =  *io  =  2  as  shown  in  Figure  3-8.  As  the  basis  functions  only  guarantee 
C°  continuity  at  u  =  1,  the  first  derivative  may  or  may  not  be  continuous  there.  Using  the 
C1  continuity  condition  (3.2),  the  first  derivative  will  be  continuous  if  and  only  if 

(<7-*4)(p§-p5)  =  (*6-*3)(p2-p§). 

Since  £4  =  t§  =  w, 

0  _  u  -  h  „p      h  -  u  _0 

*  -  —3p< +  1TT3K- 

Since  p£  =  p£  and  p°  =  P3, 

p°3  =  a3pl  +  (l-a3)p12        a3  =  ^.  (3.5) 

£7  -  13 

Note  that  (3.5)  is  the  same  as  the  deBoor  algorithm  equation  (3.4)  for  inserting  a  new  knot 
at  u  =  1,  although  the  points  and  knots  are  numbered  differently  in  this  example. 

A  similar  reasoning  yields  that  fact  that  a  knot  u  =  1  can  be  removed  a  second  time  if 
and  only  if  the  second  derivative  is  continuous,  yielding 

pJ  =  a2P2  +  (l-a2)p? 
p\  =  a3pj  +  (I  -  a3)p% 

oti  =  -   U~tl  z  =  2,3  (3.6) 

ti+p+2  ~  H 

and  a  knot  u  =  1  can  be  removed  a  third  time  if  and  only  if  the  third  derivative  is  continuous 
such  that 

P;  =  alP?  +  (l-ai)p§ 
p^  =  a2p5  +  (l-a2)Pi 

P3="3P3  +  (1-«3)P2 

a«  =      U~U  i  =  1,2,3.  (3.7) 

H+p+3  ~  H 

Note  that  there  are  no  unknowns  in  equation  (3.5),  one  unknown,  p\,  in  equation  (3.6)  and 
two  unknowns,  p\  and  p2,  in  equation  (3.7). 

For  the  knot  removal  process,  first  the  right  hand  side  of  equation  (3.5)  is  computed  and 
compared  to  p3.  If  they  are  equal  within  a  given  tolerance,  the  knot  and  p3  are  removed. 
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If  the  first  knot  removal  is  successful,  equations  (3.6)  are  solved  for  P2  and  compared: 

(l-a2)p? 


p\ 


I>2 


p\ 


Pi 


C*2 
«3P3 


1  —  as 


If  the  two  values  for  p|  are  the  same,  then  the  knot  and  control  points  P2  and  P3  are 
removed  and  control  point  P2  is  inserted. 

If  the  second  knot  removal  is  successful,  the  third  step  is  to  solve  the  first  and  third 
equations  of  (3.7)  for 


«3     _ 
Pi     = 


Pi 


Pi  -  (1  -<*i)Po 


pi 


«3P3 


1  -a3 


The  two  values  are  then  substituted  into  the  second  equation  of  (3.7).  If  the  result  is  within 
tolerance  of  P2,  then  the  knot  is  removed  and  control  points  pf,  p|  and  P3  are  replaced  by 
pf  and  pj. 


pi 


D0__1_D2. 
Pi  ~Pl  ~Pl 


Po=Po=Po=Po 


1 0  _t  1  -t  2  -t  3 


Pi 


Pi 


Pg=pi=p| 


P6=P5=P4=P3 


t4=t5=t6 


^-^-^-^o 


Figure  3-8:  Knot  Removal 

This  can  be  generalized  to  apply  to  any  number  of  removals  of  any  particular  knot.  For 
the  nth  removal,  there  will  be  a  system  of  n  equations  with  n  —  1  unknowns.  If  n  is  even, 
two  values  of  the  final  unknown  control  point  will  be  calculated  and  compared.  If  they 
are  within  tolerance,  the  knot  removal  is  successful.  If  n  is  odd,  all  new  control  points  are 
computed  and  the  final  two  are  substituted  into  the  middle  equation.  If  the  result  is  within 
the  tolerance,  the  knot  removal  is  successful.  If  the  nth  removal  is  successful,  n  control 
points  will  be  replaced  by  n  —  1  control  points. 

Knot  removal  from  a  surface  is  performed  on  the  m  +  1  rows  or  n  +  1  columns  of  control 
points.  However,  the  knot  removal  is  successful  only  if  the  knot  can  be  successfully  removed 
from  each  row  or  column.  Therefore,  the  result  must  be  checked  for  each  row  or  column 
before  any  control  points  are  removed. 
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NURBS 

The  non-uniform  rational  B-spline  (NURBS)  curve  is  represented  in  a  rational  form 

Zi=oWibiNitk(u) 


r{u) 


Zi=QmNi,k(u) 


where  W{  >  0  is  a  weighting  factor  and  Ni^u)  are  defined  over  a  non- uniform  knot  vector. 
A  NURBS  surface  patch  can  be  represented  as 


R(it,u)  = 


where  Wij  >  0  is  a  weighting  factor.  A  uniform  knot  vector  is  one  where  the  knots  are 
evenly  spaced.  If  the  knots  are  not  evenly  spaced,  the  knot  vector  is  non-uniform.  Thus, 
any  knot  vector  with  repeated  knots  is  non-uniform.  The  NURBS  representation  of  curves 
and  surfaces  allows  the  exact  representation  of  figures  such  as  circles,  conies,  quadratics, 
and  surfaces  of  revolution  with  rational  profiles  such  as  torii.  If  W{  =  1  for  all  i  or  wtj  =  1 
for  all  i,j,  this  representation  is  the  same  as  the  integral  non-uniform  B-spline  described 
previously.  In  this  thesis,  all  the  B-splines  used  are  integral,  but  the  development  could  be 
extended  to  include  NURBS. 


35 


36 


Chapter  4 

Review  of  Differential  Geometry 
Properties  of  Developable  Surfaces 

4.1  Introduction 

This  chapter  discusses  the  properties  of  ruled  and  developable  surfaces.  The  representations 
and  properties  discussed  herein  are  taken  from  classic  texts  by  Struik  [41],  Kreyszig  [26], 
Faux  and  Pratt  [11]  and  Spivak  [40]. 

4.2  Ruled  and  Developable  Surfaces 

Ruled  Surfaces 

A  ruled  surface  is  a  surface  generated  by  the  motion  of  a  straight  line  (a  generator  or  ruling) 
through  space  [41].  A  curve  that  passes  through  all  the  rulings  of  the  surface  is  called  a 
directrix.  Any  point  on  the  surface  can  be  expressed  as 

K(u,v)  =  a(u) +  v(3(u)  (4.1) 

where  a(u)  is  a  directrix  or  base  curve  of  the  ruled  surface  and  f3(u)  is  a  unit  vector  which 
gives  the  direction  of  the  ruling  at  each  point  on  the  directrix.  Alternatively,  the  surface 
can  be  represented  as  a  ruling  joining  corresponding  points  on  two  space  curves.  This  is 
represented  by 

R(u,v)  =  (1  -v)rA(u)  +vrB{u)     0<u,v<l  (4.2) 

where  rA(u)  and  r#(u)  are  directrices,  as  shown  in  Figure  4-1.  The  two  representations  are 
identical  if 

a(u)=rA(u)       and    (3(u)  =  rB{u)  -  rA(u).  (4.3) 

Throughout  the  thesis  it  is  assumed  that  the  constant  u  isoparametric  lines  correspond  to 
the  generators  of  the  developable  surface  or,  in  other  words,  the  straight  line  rulings  are  in 
the  v  direction. 

37 


directrix,  rA(u) 


generator,  ruling 


directrix,  rB  (u). 


b? . 


Figure  4-1:  A  Ruled  Surface 


Developable  Surfaces  -  A  Subset  of  Ruled  Surfaces 

Developable  surfaces  are  a  subset  of  ruled  surfaces  that  have  a  constant  tangent  plane  along 
each  generator.  Since  surface  normals  are  orthogonal  to  the  tangent  plane  and  the  tangent 
plane  along  a  generator  is  constant,  all  normal  vectors  along  a  generator  are  parallel.  This 
is  shown  in  Figure  4-2. 


tangent  plan. 


surface  normals 


Figure  4-2:  A  Developable  Surface  with  Tangent  Plane  along  a  Ruling 


Given  two  points  on  a  single  generator  that  are  selected  at  parameter  values  (uo,vi) 
and  (^0,^2)1  the  surface  is  developable  if  the  tangent  planes  at  these  points  coincide.  The 
tangent  planes  are  defined  by  the  vectors  R^  and  R^.  Using  the  representation  of  a  ruled 
surface  shown  in  (4.2),  we  obtain 


„    /           x       ,,          sdvA{uQ)  drB{u0) 

Ku(Uq,Vi)  =  (1  -Vi) —        \-Vi 


du 

Ru{uo,V2)  =  (1  -V2) ^ +V2 


du 

drB(u0) 

du 


(4.4) 
(4.5) 
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Rv{u0,vi)  =  R^uo,^)  =  rB(u0)  -rA(u0).  (4.6) 

If  the  tangent  planes  at  the  two  points  coincide,  the  vectors  Ru(uo,t>i),  Hu(uq,V2)  and 
tb(wo)  —  ta(uo)  must  be  coplanar.  Therefore,  the  triple  product  must  equal  zero: 

|Ru(wo,t;i)    Ru(u0,v2)    (rB{u0)  —  rA(u0))|  =  0.  (4.7) 

Let  Ru(uo,vi)  =  Rui  and  Ku(uq,V2)  =  R«2-  Then 

R*!  x  r„2  =  [RyulRzu2  -  RzulRyu2]i  -  [RxulRzu2  -  R-uiR^Ji  +  [R-SiR-L  -  RJiR^l*  (4-8) 

Expanding  the  first  term  of  (4.8)  by  employing  (4.4)  and  (4.5)  yields 

RyulRzu2-RzulRyu2    =    [(l-v1)i»A  +  v1i»B}[(l-v2)rzA  +  V2rB] 

-[(1  -  vx)vA  +  Wlrf,][(l  -  v2)v\  +  v2ryB] 
=     (1  -  i/xXl  -  W2)fjfi  +  «i(l  -  v2)r^  +  W2(l  -  ^i)r^r2B  +  vxv2ryBrzB  - 

[(1  -  Vl)(l  -  v2)t\va  +  «i(1  -  t^rf,  +  t*(l  -  Wl)r*rJ  +  ^2rysr2B] 
=     (vi  -  v\v2  -v2  +  v\v2)rzAryB  +  {v2  -  v\v2  -  v\  +  viV2)ryAtB 
=     {v2-vi)(ryArzB-rAryB). 


Similarly, 


and 


Therefore, 


R^R22  -  R^R^  =  (v2  -  wi)(r^r^  -  r^r£) 


RS1RS2  -  R?i^2  =  («2  -  ^(^4  -  *M). 


Rul  X  R^  =  (V2  -  vi)(ta  x  r5) 

and  thus  (4.7)  can  be  reduced  to 

(rB(u)  -  r^(w))  x  rA{u)  ■  rB{u)  =  0  (4.9) 

for  any  u  as  long  as  v\  7^  i^-  Conversely,  any  surface  that  satisfies  (4.9)  must  satisfy  (4.7) 
and,  therefore,  the  tangent  planes  must  be  the  same  and  the  surface  must  be  developable. 
Thus,  a  surface  is  developable  if  and  only  if  (4.9)  is  satisfied.  Substituting  (4.3)  into  (4.9) 
and  using  the  fact  that  r.4  x  r^  =  0  yields  the  equivalent  condition 

ax  (3-0  =  0.  (4.10) 

Developable  Surfaces  -  Envelopes  of  Families  of  Planes 

An  alternative  representation  of  a  developable  surface  is  as  an  envelope  of  a  family  of  planes. 
The  concept  of  an  envelope  will  be  described  here  for  a  family  of  curves,  then  extended  to 
planes.  First,  we  represent  a  curve  as  r(u,  a)  where  u  is  the  curve  parameter  and  a  is  the 
family  parameter.  In  other  words,  each  curve  in  the  family  is  v(u,  an).  There  may  be  a 
curve  re(a)  that  is  tangent  to  each  curve  in  the  family,  as  shown  in  Figure  4-3.  This  curve 
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would  be  the  envelope  of  that  family  of  curves  [11]. 


r  -  r. (a) 


r  -  rCu.ap 


r  -  r(u,a,) 


Figure  4-3:  An  Envelope  of  a  Family  of  Curves 

Similarly,  a  surface  can  be  the  envelope  of  a  family  of  surfaces.  Let  a  one-parameter 
family  of  surfaces  be  represented  by  G(x,a),  where  a  is  the  family  parameter,  and  assume 
that  the  surfaces  are  twice  continuously  differentiable.  Also  assume  that  consecutive  surfaces 
intersect  one  another.  The  intersection  is 

G(x,a)=0       G(x,a  +  /i)=0 

where  \h\  is  sufficiently  small.  This  may  also  be  represented  by 

G{x,a  +  h)-G(x,a)  =Q 
h 

The  limit  of  (4.11)  as  h  tends  to  zero  is  dG(x,  a) /da  —  0.  The  set  of  points  represented  by 

G(x,Q)  =  o     ^)=0 

oa 

is  the  characteristic  of  the  surface.  If  the  characteristics  of  all  the  possible  surfaces  in  the 
family  form  a  surface,  that  surface  is  called  the  envelope.  If  the  envelope  exists,  then  at 
every  point  of  the  characteristic  of  a  surface  in  the  family,  the  tangent  planes  of  the  surface 
and  of  the  envelope  coincide. 

A  characteristic  point  is  determined  by  the  intersection  of  the  characteristic  with  another 
surface  in  the  family.  These  points  of  intersection  are  determined  by 

G(x,a)=0        dG^a)  =Q       G(x,a  +  A;)=0. 
oa 

where  \k\  is  sufficiently  small.  Using  the  Taylor  series  expansion 

r%      rM       s      ,dG{x,a)       k2d2G(x,a  +  9k)        n       _      , 
G(x,  a  +  k)  =  G(x,  a)  +  k — ^-!-  +  — ^ "'       0  <  »  <  1, 

40 


the  third  equation  can  be  rewritten  as 

d2G(x,a  +  9k) 
da2 


=  0. 


Then,  taking  the  limit  as  k  approaches  zero,  the  characteristic  points  of  the  surface  are 


G(x,  a)  =  0        — - =  0 


d2G(x,a) 
da2 


=  0. 


If  the  characteristic  points  exist  and  if  they  form  a  curve  as  a  varies,  the  curve  lies  on  the 
envelope  of  the  family  and  is  called  the  edge  of  regression  of  the  envelope.  If  the  envelope 
and  the  edge  of  regression  exist,  then,  at  every  point  of  the  edge  of  regression,  the  tangents 
of  the  edge  of  regression  and  the  corresponding  characteristics  coincide.  The  characteristics 
then  form  two  sheets  corresponding  to  the  tangents  of  the  edge  of  regression  in  the  positive 
and  negative  directions.  These  sheets  form  a  cusp  at  the  edge  of  regression,  as  shown  in 
Figure  4-4. 


Second    She 


Figure  4-4:  Edge  of  Regression 
Let  us  now  consider  a  family  of  planes  and  their  derivatives  which  can  be  represented 


by 


C(x,a) 

dG 

da 

d2G 

da2 


a(a)  •  x  +  c(a)  =  0 
da 


dc_ 
da  da 


d2a 
da2 


x  + 


d2c 
da2 


0 


=  0 


where  a(a)  is  the  unit  normal  vector  to  the  corresponding  plane.  When  all  planes  are  paral- 
lel, no  intersection  exists  and  no  characteristic  curves  are  formed.  When  all  planes  intersect 
in  a  single  line  forming  a  pencil,  the  envelope  is  a  line  instead  of  a  surface.  Therefore,  these 
two  cases  are  excluded  [26].  When  the  envelope  exists,  it  is  formed  of  the  characteristics  of 
the  family  of  planes.  Since  planes  intersect  in  straight  lines,  the  envelope  is  a  ruled  surface. 
Furthermore,  the  tangent  plane  along  a  characteristic  is  constant  since  this  is  an  envelope 
of  planes,  and  hence  the  surface  is  developable. 
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When  the  vectors  a(a)  are  all  parallel  to  a  single  plane  P,  then  all  the  planes  in  the 
family  are  parallel  to  a  vector  normal  to  the  plane  P.  Thus,  the  envelope  of  the  family  of 
planes  forms  a  cylinder  (except  when  they  form  a  pencil  which  has  been  excluded).  When 
all  the  characteristic  points  coincide,  all  the  characteristics  intersect  in  a  single  point  and 
the  surface  is  a  cone.  When  the  characteristic  points  form  a  curve,  the  curve  is  the  edge 
of  regression  and  the  surface  is  the  tangent  surface  to  the  curve.  Therefore,  all  developable 
surfaces  can  be  represented  as  planes,  cylinders,  cones  and  tangent  surfaces,  or  combinations 
thereof  [41]. 

This  description  of  developable  surfaces  leads  to  the  following  duality  between  space 
curves  and  developable  surfaces  [41,  p. 72]: 


Curves 

2  points  determine  a  line. 

3  points  determine  a  plane. 
2  consecutive  points  on  a 

curve  determine  a  tangent  line. 

The  curve  is  the  envelope 
of  tangents. 


Developable  Surfaces 

2  planes  determine  a  line. 

3  planes  determine  a  point. 
2  consecutive  planes  of  the 

family  of  planes  determine 
a  characteristic  line. 
The  developable  surface  is 

generated  by  characteristic  lines. 


Gaussian  Curvature 

Using  equation  (4.1) 

R(u,  v)  =  a(u)  +  v/3(u), 

the  first  and  second  derivatives  of  a  surface  can  be  represented  as 

R^tt,?;)  =  a(u)+v0(u) 

Ry(u,v)  =  (3(u) 

RuU(u,v)  =  ct(u)  +  v(3(u) 

Ruw(u,u)  =  /3(u)      and 

RyV(u,v)  =  0. 

Also  recall  from  (2.9)  that 

EG-F2  =  |Ru  x  R„|2  >0 

for  a  regular  surface.  The  surface  normal  is  given  by  (2.7) 

R-u  x  R-v 


N    = 


|Ru  x  R^l 
{a{u)+v0{u))  x  0(u) 

VEG-F2 
a{u)  x  (3{u)  +  v[/3(u)  x  (3(u)) 

y/EG  -  F2 
a(u)  x  /3(u) 

"Teg  -  f2  ' 


(4.12) 

(4.13) 
(4.14) 
(4.15) 
(4.16) 


(4-17) 
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Recall  from  (2.23)  that  Gaussian  curvature  is  represented  as 

LN  -M2 
EG-F2' 

where,  from  equations  (2.13)  and  (2.12) 

N    =    KyV  ■  N     and 
M    =    R^-N. 

Thus,  from  (4.16),  N  =  0  and  the  Gaussian  curvature  can  be  rewritten  as 

M2 


K 


EG-F2 
(Ru?  •  N)2 
EG-F2 

(/3(u)  ■  a(u)  x  (3(u))2 
(EG  -  F2)2 

(QX/3-/3)2 
(EG-F2)2' 


(4.18) 


Thus,  from  condition  (4.10),  zero  Gaussian  curvature  is  a  sufficient  and  necessary  condition 
for  a  surface  to  be  developable. 

Since  the  Gaussian  curvature  of  a  developable  surface  is  zero  everywhere,  the  maximum 
and  minimum  principal  curvatures  (2.21)  and  (2.22)  of  a  developable  surface  can  be  written 
as 

Kmax  =  -"  ~r  I-"  l>      Kmin  =  •"        |-"r 

These  principal  curvatures  reduce  to 

Kmax  =  2H,      Kmin  =  0      when  H  >  0,  (4.19) 

Kmax  =  0,      Kmin  =  0      when  H  =  0     and  (4.20) 

Kmax  =  0,    Kmin  =  2#    when  H  <  0.  (4-21) 

It  is  clear  from  equations  (4.19)  to  (4.21)  that  at  least  one  of  the  principal  curvatures  is 
always  zero  on  a  developable  surface,  which  agrees  with  the  fact  that  the  Gaussian  curvature 
is  zero  everywhere  (see  equation  (4.18)).  Kmax  and  Kmin  from  (4.19)  and  (4.21)  respectively 
are  each  termed  the  nonzero  principal  curvature,  k*,  where 


k*  =  2H.  (4.22) 


Mapping 


Isometric  mapping  is  defined  as  a  mapping  that  preserves  arc  lengths  between  the  mapped 
surfaces.  Two  surfaces  that  can  be  isometrically  mapped  to  one  another  are  called  isomet- 
rics. Since  arc  lengths  are  preserved,  the  first  fundamental  forms  of  isometrics  are  the  same. 
Since  the  Gaussian  curvature  of  a  surface  and  the  geodesic  curvature  of  a  curve  on  a  sur- 
face depend  only  upon  the  coefficients  of  the  first  fundamental  form,  corresponding  points 
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on  two  isometric  surfaces  have  the  same  Gaussian  curvature  and  corresponding  curves  on 
isometric  surfaces  have  the  same  geodesic  curvature  at  corresponding  points  [26]. 

A  conformal  mapping  is  one  in  which  the  angle  between  every  pair  of  intersecting  arcs 
is  the  same  on  both  surfaces.  This  occurs  if  the  coefficients  of  the  first  fundamental  form  on 
the  original  surface  are  proportional  to  those  on  the  inverse  surface  [26].  Therefore,  every 
isometric  surface  is  conformal. 

Since  a  plane  has  zero  Gaussian  curvature,  only  a  surface  with  zero  Gaussian  curvature 
can  be  mapped  to  it.  Therefore,  only  a  developable  surface  can  be  isometrically  mapped  to 
a  plane  [26,  p.  189]. 

Since  the  tangent  planes  along  any  ruling  are  constant,  developable  surfaces  can  be  un- 
rolled, or  developed,  isometrically  into  a  plane  without  stretching  or  tearing  by  successively 
placing  each  ruling  on  the  plane.  One  can  think  of  a  developable  surface  as  a  paper  with 
folds.  As  the  number  of  folds  approaches  infinity,  the  surface  becomes  smooth.  The  surface 
can  then  be  unfolded  again  into  a  plane.  This  is  shown  in  Figure  4-5. 


Figure  4-5:  A  Surface  with  Folds  and  a  Developable  Surface 

A  geodesic  represents  the  curve  with  the  shortest  length  between  two  points  on  a  sur- 
face. Therefore,  geodesies  on  a  plane  are  straight  lines.  Since  developable  surfaces  map 
isometrically  to  planes,  geodesies  on  developable  surfaces  map  to  straight  lines  on  a  plane 
and  straight  lines  on  a  plane  map  to  geodesies  on  a  developable  surface. 

Developable  Surface  Properties 

The  following  statements  are  equivalent  necessary  and  sufficient  conditions  for  a  surface  to 
be  developable.  Each  has  been  described  above.  Thorough  proofs  can  be  found  in  Struik  [41] 
and  Kreyszig  [26]. 

1.  Developable  surfaces  possess  the  same  tangent  plane  at  all  points  of  the  same  gener- 
ator. 

2.  The  normal  vectors  on  a  developable  surface  along  a  ruling  are  parallel. 

3.  A  developable  surface  is  the  envelope  of  a  one-parameter  family  of  planes. 

4.  Gaussian  curvature  is  zero. 
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5.  The  mapping  of  a  developable  surface  onto  a  plane  is  isometric. 

6.  Geodesies  on  a  developable  surface  map  to  straight  lines  on  a  plane,  and  straight  lines 
on  a  plane  map  into  geodesies  on  a  developable  surface. 
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Chapter  5 

Design  With  B-Spline  Developable 
Surfaces 

5.1     Introduction  and  Literature  Review 

In  the  current  literature,  there  are  three  main  approaches  to  representing  developable  sur- 
faces in  Bezier  form.  These  approaches  are  described  below.  This  chapter  explores  user- 
friendly  methods  for  the  design  of  developable  surfaces.  Section  5.2  describes  the  design  of 
developable  B-spline  strip  surfaces  that  lie  between  parallel  planes.  Section  5.3  describes 
the  design  of  developable  B-spline  surfaces  with  general  three-dimensional  directrices.  The 
development  of  a  developable  surface  onto  a  plane  is  described  in  Section  5.4.  This  is 
an  efficient  implementation  of  the  method  in  the  current  literature,  also  described  below. 
Section  5.5  provides  some  engineering  examples. 

5.1.1     Developable  Bezier  Strips 

Aumann  [1]  provides  the  basis  for  a  user  friendly  and  computationally  inexpensive  method 
of  designing  developable  surfaces.  He  represents  developable  surfaces  in  terms  of  two  Bezier 
curves  (directrices)  and  rulings  between  pairs  of  points  from  each  curve,  as  discussed  in 
Section  4.2,  but  restricts  the  two  directrices  (ta(u),tb(u))  to  lie  on  parallel  planes  so  that 
the  tangent  vectors  ta(u)  and  tb(u)  are  parallel,  i.e.  r#(u)  =  p{u)ta{u),  where  p(u)  is  a 
linear  function  of  u.  This  condition  automatically  satisfies  the  condition  (4.9). 

His  design  philosophy  is  to  consider  r^(u)  as  a  design  curve,  which  means  its  degree,  knot 
vector  and  control  points  are  specified.  In  addition,  the  two  end  points  r#(0)  and  r#(l)  of 
the  resulting  directrix  are  provided  as  boundary  conditions.  The  knots  and  internal  control 
points  of  the  second  directrix  rg(u)  are  then  computed  using  a  simple  system  of  linear 
equations.  He  investigates  the  case  when  the  design  curve  is  a  cubic  Bezier  curve  and  p(u) 
is  a  linear  function  of  u.  Hence,  the  resulting  directrix  is  quartic.  Figure  5-1  illustrates  the 
terminology. 

Frey  and  Bindschadler  [14]  extend  the  work  of  Aumann  [1]  by  generalizing  the  degree  of 
the  directrices  and  by  considering  the  case  where  p(u)  is  quadratic.  Aumann  [2]  further  ex- 
tends his  work  so  r^(u)  and  r#(u)  have  the  same  degree.  All  three  papers  discuss  regularity 
conditions  of  the  developable  surfaces. 
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design  curve.  rA(u) 


generator,  ruling 


resulting  directrix,  rg  (u) 


D21 


Figure  5-1:  A  Developable  Surface 

5.1.2  Duality  Between  Points  and  Planes 

The  second  approach,  introduced  by  Bodduluri  and  Ravani  [3],  is  based  on  the  represen- 
tation of  developable  surfaces  as  an  envelope  of  a  one-parameter  set  of  planes  as  described 
in  Section  4.2.  They  express  the  surfaces  in  terms  of  plane  geometry  using  the  concept  of 
duality  between  points  and  planes  in  3D  projective  space.  They  interpret  the  developable 
surface  as  the  set  of  its  tangent  planes  represented  in  the  dual  form.  Their  method  is  math- 
ematically elegant;  however,  it  is  somewhat  less  user  friendly  and  requires  a  conversion  to 
a  standard  tensor  product  form  [37]. 

Pottmann  and  Farin  [37]  generalize  the  results  of  Bodduluri  and  Ravani  [3]  by  providing 
algorithms  for  converting  the  dual  representation  of  developables  to  standard  tensor  product 
form.  Hoschek  and  Pottmann  [23]  derive  algorithms  for  interpolation  or  approximation 
with  developable  B-spline  surfaces.  Given  a  developable  surface  patch  ft,  they  pick  out 
an  appropriate  number  of  tangent  planes  Y*  of  a  developable  surface  patch  ft,  and  either 
interpolate  or  approximate  this  data  with  a  developable  B-spline  surface. 

5.1.3  Nonlinear  Representation 

The  third  approach  is  due  to  Lang  and  Roschel  [29].  They  obtain  conditions  for  the  control 
nets  and  weights  of  rational  Bezier  surfaces  to  be  developable,  which  leads  to  a  complicated 
nonlinear  system  of  conditions.  Although  their  work  provides  a  useful  criterion  for  checking 
the  developability  of  a  given  ruled  surface  patch,  it  is  difficult  to  design  a  developable  surface 
with  this  method. 


5.1.4     Development 

Faux  and  Pratt  [11]  provided  a  method  for  the  development  of  a  curve  on  a  developable 
surface  onto  a  plane  based  on  the  fact  that  a  developable  surface  maps  isometrically  onto  a 
plane.  This  method  is  also  described  by  Gurunathan  and  Dhande  [16].  Another  method  is 
described  by  Clements  and  Leon  [7],  in  which  a  geodesic  on  the  surface  is  mapped  to  the 
plane  and  points  on  the  plane  are  determined  in  reference  to  this  geodesic.  This  method, 
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however,  requires  the  mapping  of  the  geodesic  and  then  the  calculation  of  the  location  of 
each  edge,  and  is  therefore  slower  than  the  other  method  [11]. 


5.2     B-Spline  Developable  Surfaces  with  Directrices  in  Paral- 
lel Planes 

This  section  further  extends  Aumann's  work  [1]  such  that  n  Bezier  developable  surfaces, 
whose  two  directrices  lie  in  parallel  planes,  are  strung  together  by  joining  them  along  their 
end  rulings.  In  this  manner,  they  can  be  represented  as  a  single  B-spline  surface  with  C2 
continuity. 

From  Section  4.2,  a  surface  is  developable  if  and  only  if  the  tangent  plane  is  constant 
along  a  generator.  If  the  two  directrices  ta(u)  and  tb(u)  are  restricted  to  lie  in  parallel 
planes,  then  the  tangent  vectors  to  r^(u)  and  rs(n)  must  be  parallel  at  every  u  for  the 
surface  to  be  developable.  Therefore,  condition  (4.9)  is  equivalent  to 

rB(u)  =  p(u)rA{u)  (5.1) 

where  p(u)  is  a  scalar  function  of  u.  In  this  thesis,  p{u)  is  selected  to  be  linear,  expressed 
as 

p{u)  =po(l  -w)  +  Piu. 

The  design  philosophy  is  to  consider  r^(u)  as  a  design  curve  [1,  14],  which  means  its 
degree,  knot  vector  and  control  points  are  specified.  Also,  the  two  end  points  r5(0)  and 
r#(l)  are  provided  as  boundary  conditions.  The  knots  and  internal  control  points  of  the 
second  directrix,  r#(u),  are  then  computed.  Since  multiplication  of  B-spline  functions  is 
complex  and  computationally  expensive,  the  Bezier  segments  are  first  extracted  from  the 
B-spline  curve  by  knot  insertion  using  Boehm's  algorithm  [22]  described  in  Section  3.3. 

Figure  5-2 (a)  shows  the  design  curve  r^(u)  which  is  a  planar  cubic  B-spline  curve  with 
knot  vector  T  =  (0  0  0  0  ±  §  1  1  1  1)  and  control  points  (0,  2,  0),  (0.8,  0.5,  0),  (1.6,  1.5, 
0),  (2.4,  1.5,  0),  (2.6,  -0.5,  0)  and  (4,  0,  0).  Figure  5-2(b)  shows  three  Bezier  segments 
after  knot  insertion  with  the  two  endpoints  of  tb(u),  which  are  supplied  by  the  user  to  be 
rs(0)  =  (  —  1, 3,  —3)  and  rg(l)  =  (2, 1.3,  —3).  Note  that  ta(u)  is  in  the  z  =  0  plane  and  the 
tb  endpoints  are  in  the  z  =  3  plane,  parallel  to  the  plane  of  r^. 

Let  ta{u)  be  a  degree  n  curve  with  m  Bezier  segments.  Since  p(u)  is  linear,  the  degree 
of  tb(u)  is  n  -I-  1.  We  denote  the  jth.  Bezier  segments  as 

n 

^(w)     =     5Z  ajn+i#tin(u)       j  =  0, . . . ,  m  -  1  and 

n+l 

rflfa)     =     J2h](n+i)+iBi,n+i{u)      j  =  0,...,m-l.  (5.2) 

The  control  points  a  and  b  are  two-dimensional  vectors  since  they  have  been  restricted  to 
lie  in  parallel  z  planes.  Substituting  (5.2)  into  equation  (5.1)  yields 

n 

(n  +  1)  5^(bj(n+i)+i+1     -    bj{n+l)+i)Bi<n(u) 

i=Q 
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(a)  Design  Curve  r.(u)  with  Control  Polygon  and  Knots 


(b)  Design  Curve  after  Knot  Insertion  with  Bezier  Segments 
and  Endpomts  of  rB(u)  Shown 


(c)  Design  Curve  and  Resulting  Directrix  before  Knot  Removal  (d)  Final  Developable  B-Spline  Strip  Surface 

Figure  5-2:  Design  of  Developable  B-Spline  Strip  Surface 


n-i 


=     p(u)n^2(Sijn+i+i  -  SLjn+i)Bi7n-i{u) 

n-1 

=    n^2(p0(l  -  u)  +  piu)(ajn+i+i  -  aijn+l)Btjn-i(u)    (5.3) 


i=0 


for  j  =  0, . . . ,  m  —  1.  Using  the  definition  of  a  Bernstein  polynomial  (3.1), 


(l-u)fliin_i(«)     =     (1-u) 


1)!  (i-.r-^-v 


i!(n-l  -*)! 
n  -  1      (n)' 


n     i\(n  —  i)\ 

n  —  1  _     .   . 
■BinW 


(l-u)n-lul 


n 


(5.4) 
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and 


n     {i  +  l)!(n  -  {i  +  l)!v  ; 

^i+i.nM-  (5-5) 


n 
Substituting  (5.4)  and  (5.5)  into  (5.3)  yields 


n 

(n  +  1)  5^(bj(n+l)+»+l      ~      hj(n+l)+i)Bi,n(u) 
i=0 

n-1 

=      n  S  f^1  ~  u)-Bi,n-l(w)  +  PluBi,n-l{u)]  (ajn+j+i  -  BLjn+i) 
i=Q 
n-1 

=      ^2  t^°(n  ~  i)Bi,n(u)  +  Pl(*  +  l)^i+l,n(w)]  (ajn+i+i  -  ajn+i) 


i=0 

pon(ajn+i  -  ajn)Bo,n(w) 

n-1 
+  5Z  ^°(n  ~  0(ai"+i+l  ~  ajn+i)  +  Pii(aj„+j  -  aLjn+i-i)]  Bitn(u) 


t=l 


+pin(&jn+n  -  &jn+n-i)Bnfn(u). 

Since  Bernstein  basis  polynomials  are  linearly  independent,  each  set  of  coefficients  must  be 
equal,  yielding  the  following  system  of  n  —  1  equations  [14]  for  the  jth  Bezier  segment 

np^Ajn+i  -  &jn)     =     (n  +  l)(bj(n+1)+1  -bj(n+1))       (5.6) 
(n  -  i)pj(aj„+j+i  -  a.jn+i)  +  ip{(&jn+i  -  ajn+i-i)     =    (n  +  l)(bJ,(n+1)+j+1  -  bj(n+1)+t) 

n^i(a0'+i)n  _  a(j+1)n_1)     =     (n  +  l)(b(j+1)(n+1)  -  by+i)(n+i)-i)i 

where  i  =  l,...,n  —  1. 

Now  consider  joining  two  consecutive  Bezier  patches  with  C2  continuity.  First,  C1 
continuity  requires  that  the  first  derivatives  along  the  connecting  edges  of  both  patches 
must  be  equal.  Since  the  connecting  edge  is  a  ruling  and  is  therefore  linear,  this  requirement 
reduces  to 

r^(l)     =    r^+1(0)  (5.7) 

*jj(l)     =    ^+1(0).  (5-8) 

From  equation  (5.1)  we  have 

ri(l)    =    ^'(l)ri(l) 

ii+1(0)     =    p>"+1(0)ri+1(0).  (5.9) 

Substituting  equations  (5.9)  into  equation  (5.8)  and  using  equation  (5.7)  yields 

pi=V0+1-  (5-10) 
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C2  continuity  requires  that  the  second  derivatives  along  the  connecting  edges  must  be 


equal.  Again,  the  connecting  edge  is  linear,  which  reduces  this  requirement  to 

r^(l)    =    i^+1(0) 
*i(l)     =    ^+1(0)- 
Taking  the  derivative  of  equation  (5.1)  yields 

rfl(u)  =  p(u)rA{u)  +  p(u)tA(u) 

where 

p(u)  =pi  -po. 

Evaluating  equation  (5.13)  at  the  connecting  edge  yields 

^(0)     =    ^o+1^+1(0)  +  (pi+1-^0+1)r^1(0). 
Substituting  equations  (5.14)  into  equation  (5.12)  and  using  equation  (5.11)  yields 


(5.11) 
(5.12) 


(5.13) 


f+i 


2p>0+l     =    pJ  +  />{+1. 


(5.14) 


(5.15) 


Each  Bezier  patch  has  n  +  1  control  points  along  rA(u)  and  n  +  2  control  points  along 
r#(u),  where  n  is  the  degree  of  the  original  B-spline  design  curve.  The  system  of  equations 
for  a  single  patch  derived  in  equation  (5.6)  can  be  solved  for  the  coefficients  of  the  second 
directrix 

bj(n+l)+l      =      bj(n+1)  +  — — -po(aj'n+l  -  ajn) 

,  ,  n  -  i  j 

Dj(n+l)+i+l      —      Bj(n+l)+i  ^ ~ ~7Po\aJn+i+l  ~  &jn+i) 

i 

H — rp{{a-jn+i  -  ajn+j-i),  »  =  l,...,n  — 1 

n+  1 

n 

b(j  +  l)(n+l)      =      b(j  +  l)(n+l)-l  +  ^7YPl(a(j  +  l)n  -a(j  +  l)n-l)  (5-16) 

where  j  =  0, . . . ,  m  —  1  for  a  system  of  m  patches.  Adding  equations  (5.16)  for  m  patches 
yields  a  vector  equation  that  can  be  solved  for  p§  and  p™~x  in  terms  of  the  first  and  last 
points  of  the  second  directrix: 


m—  1, 


bm(n+l)  _  b0      =     PovO  +  P™      vl 


(5.17) 


where 


vo 


vi     = 


1_ 

m(n  +  1)  V^o 


m— 1 


/run 

(Ea*  +  E  »*" 
i=l  y 


ao 


and 


&mn 


1 

m(n  +  1) 


m— 1 


E^+E 


u=0 


i=l 
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Combining  equations  (5.10)  and  (5.15),  we  can  express  pQ  and  p\  in  terms  of  the  two 
end  coefficients  p[j  and  p™-1  as  follows: 

Therefore,  given  a  design  directrix,  ta(u),  and  the  endpoints  of  the  second  directrix,  1:5(0) 
and  r^(l),  all  the  interior  control  points  of  the  second  directrix  can  be  obtained  succes- 
sively using  equations  (5.16).  After  the  second  directrix  is  calculated,  the  first  directrix  is 
degree  elevated  so  the  two  curves  may  be  combined  into  a  B-spline  strip.  Any  knots  that 
may  be  removed  from  the  surface  while  leaving  the  curve  geometrically  and  parametrically 
unchanged  are  removed  using  the  knot  removal  technique  [42]  described  in  Section  3.2. 

Figure  5-2(c)  shows  the  computed  tb(u)  together  with  r^u)  in  Bezier  form  and  the  two 
end  rulings.  Figure  5-2 (d)  shows  the  final  B-spline  developable  surface  after  knot  removal. 

Developable  surfaces  become  irregular  when  |Ru  x  RJ  vanishes,  which  leads  to  the 
regularity  condition  for  a  strip  surface  [14] 

p{u)  >  0,  0  <  u  <  1. 

From  a  user's  point  of  view,  once  he/she  determines  one  of  the  endpoints,  it  is  convenient 
to  know  the  range  that  the  other  end  point  can  take.  From  (5.18)  it  can  be  seen  that 
if  both  Pq  and  p™~  are  positive,  then  all  p  coefficients  are  positive  and  the  surface  is 
regular.  The  coefficients  p  can  be  determined  by  solving  the  vector  equation  (5.17).  The 
resulting  equations  for  p[j  and  Pi"-1  must  be  positive  for  the  surface  to  be  regular,  yielding 
the  following  restrictions: 

o  _  W»(„+1)  -  jgj  -  (%(„+1)  -  gK 


v%v\  -  v^v0 


and 


,  =  (^o-gM-Wy+o-WS  >  0 

VqV\  —  Vf  Vq 

where  the  superscripts  denote  the  x  and  y  components  of  the  vectors.  Given  one  endpoint, 
a  regular  surface  may  be  formed  by  selecting  the  second  endpoint  to  lie  within  the  range  of 
values  determined  by  equations  (5.19)  and  (5.20). 

It  is  also  interesting  to  observe  [14]  that  when  the  two  vectors  joining  two  end  points  of 
each  directrix  are  parallel  and  their  magnitudes  are  multiples  of  a  constant,  the  resulting 
developable  surface  is  a  generalized  cone.  When  the  constant  is  equal  to  one,  the  surface  is 
a  cylindrical  surface.  In  other  words,  if  p(u)  is  a  constant  not  equal  to  one,  the  surface  is  a 
generalized  cone,  and  if  p(u)  =  1,  the  surface  is  a  cylinder. 

Algorithm 

The  following  steps  describe  the  algorithm  for  the  solution  method. 

1.  The  user  provides  the  design  curve  r^(tt)  and  the  endpoint  i"b(0).  The  program 
provides  a  range  of  values  for  r#(l)  that  ensure  a  regular  surface,  and  the  user  supplies 
a  specific  value  for  rfl(l). 
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2.  Bezier  segments  are  extracted  from  r^(u)  by  knot  insertion. 

3.  The  system  of  equations  (5.16)  is  solved  for  the  internal  control  points  of  r#(u). 

4.  Once  the  resultant  directrix  r#(it)  is  computed,  the  developable  Bezier  surface  is 
obtained  by  degree  raising  r^(w)  and  substituting  in  equation  (4.2). 

5.  Knots  that  can  be  removed  while  leaving  the  surface  unchanged  are  removed  using 
the  method  described  in  Section  3.2. 

5.3     Developable  Surfaces  with  3D  Directrices 

5.3.1     Developable  (n,  1)  Bezier  Surfaces 

In  Section  5.2,  we  restricted  the  two  directrices  to  lie  on  parallel  planes.  In  such  case,  r^(u) 
and  tb(u)  must  be  parallel  to  satisfy  condition  (4.9).  In  this  section,  the  directrices  are  not 
constrained  to  lie  on  parallel  planes,  so  equation  (5.1),  tb  =  p(u)ta,  no  longer  applies. 

The  design  philosophy  here  is  the  same  as  that  of  the  previous  section.  Namely,  r^(u) 
is  designated  the  design  curve  with  its  degree,  knot  vector  and  control  points  specified.  The 
two  end  points  r#(0)  and  rs(l)  are  provided  as  boundary  conditions  and  the  knots  and 
internal  control  points  of  the  second  directrix,  tb{u)  are  computed. 

We  begin  this  development  by  assuming  the  design  curve  ta(u)  is  a  Bezier  curve  of 
degree  n, 

n 

rA{u)  =  5^aj.Bjin(u),        0  <  u  <  1, 

i=0 

where  az  =  (af,a^,a^)T  are  control  points.  The  resultant  curve  r#(u)  is  also  a  Bezier  curve 
of  degree  n.  In  Section  5.3.2  we  investigate  the  case  when  the  curves  r^(w)  and  tb{u)  are 
B-spline  curves.  If  bj  =  (bf,bf,b*)T,  i  =  0, ...  ,n,  are  the  control  points  of  the  degree  n 
Bezier  curve  tb{u),  and  a*,  b^  and  Cj  are 

aj  =  (af,af,a|)r  =  al+i -aj  z  =  0,  ...,n-l 
b,  =  (6f,6f,6t2)T  =  bl+1-bl  *  =  0,...,n-l 
Ct     =     (cf,cf,ctz)T  =  bt -aj  z  =  0,  ...,n, 

then  equation  (4.9)  can  be  expressed  as 

(jT^niu))  x  (n^a^-itw))  •  (n^2btBltn_l{u)\  =0.  (5.21) 

Since  equation  (5.21)  is  a  scalar  polynomial  equation  of  degree  3n  —  2,  we  can  rewrite  it 
(upon  division  by  n2)  in  the  form 

3n-2 

£  EkBk>3n-2(u)  =  0  (5.22) 

Jfc=0 
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where 

min(2n—l,k)  i2n— lwn- 1\ 

Ek=     y:     (p^t-i-prH-i+PiH-i)  isn-t*         (5-23) 

i=max(0,k— n+l)  \    k    ' 

and 

min(n,i)  (n\  (n~^\ 

Px  =  V  (cya*     -czay    rJ     l~J 

*t    —  Z^  Vcjai-j        c]ai-j)     (2n-l\ 

j=max(Q,i—n+l)  v     z     / 

77iin(n,x)  /n\  /n— 1\ 

*t  Z^  ^Cj"t-J        C3ai-3>     /2n-l\ 

j=max(0,i— n+l)  V     i     ' 

mm(n,t)  /n\  m— 1\ 

P2  =  V  (cxay     -cyax    )J     l~J 

^*  ZL  vCJai-J        C]ai-3>     /2n-l\ 

j=max(0,i—n+l)  \     i     ' 

using  the  algorithms  for  multiplication  of  polynomials  in  Bernstein  form  described  by 
Farouki  and  Rajan  [10].  Since  equation  (5.22)  is  valid  for  all  u  and  the  Bernstein  basis 
functions  are  linearly  independent,  Ek  must  be  zero  for  all  k,  0  <  k  <  3n  —  2.  This  leads  us 
to  a  system  of  3n  —  1  polynomial  equations  with  3(n  —  1)  unknowns,  i.e.  bj  (i  =  1, . . . ,  n  —  1) 
where  the  two  boundary  points,  bo  and  bn,  are  given.  Therefore,  there  are  two  more  equa- 
tions than  unknowns  and  the  system  is  overdetermined.  The  system  of  equations  was  solved 
using  a  Numerical  Algorithms  Group  program  [35]  for  minimization. 

Solution  Method  Using  Unconstrained  Optimization 

One  way  to  solve  the  overdetermined  system  is  to  seek  a  least  squares  solution  to  E(x)  =  0 
where 

X     =      (blbl,bl,...X-l,byn-lX-l)T 
and 

E(x)     =     (£0(x),£i(x),...,£3n-2(x))T. 

This  can  be  obtained  by  minimizing  the  function 


1  11  3n— 2 

F(x)  =  -  (£02(x)  +  £?(x)  +  . . .  +  i&_2(x))  =  2ETWE(X)  =  2  £  ^(x)'      (5"24) 

The  problem  can  then  be  reduced  to  searching  for  zeros  of  the  gradient  vector  field  VF(x), 
i.e. 

VF(x)  =  0.  (5.25) 

The  Newton- Raphson  method  can  be  used  to  iteratively  solve  the  nonlinear  system  (5.25). 
The  Taylor  series  expansion  about  the  current  point  is 

VF(x  +  5x)  «  g(x)  +  [H(x)]<5x  (5.26) 
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where  g(x)  and  [H(x)]  are  the  gradient  vector  and  the  Hessian  matrix  of  the  objective 
function  F(x).  Using  (5.26),  the  step  direction  Sx  is  found  from 

[H(x)]<Jx  =  -g(x).  (5.27) 

The  value  for  the  subsequent  iteration,  x  is 

x  =  x  +  Sx. 

Using  the  special  properties  of  a  least  squares  function  such  as  equation  (5.24),  the  deriva- 
tives of  .F(x)  are  represented  as 

g(x)  =  VF(x)  =  [A(x)]E(x) 

and 

3n-2 

[H(x)]  =  V2F(x)  =  [A(x)][A(x)]T  +  £  ^(x)V2^(x) 

where 

[A(x)]  =  [VEo,V£a,...,V£3n-2] 

is  the  (3n  —  3)  x  (3n  —  1)  Jacobian  matrix,  the  columns  of  which  are  the  first  derivative 

BE 

vectors  VjE^  of  the  components  of  E  ([Ay]  =  tjt?-)- 

In  any  minimization  process,  the  initial  approximation  for  x  greatly  affects  the  resulting 
solution.  Here  we  use  a  linear  interpolation  to  obtain  the  initial  values  for  b^,  where 

i                       tx  —  i 
bz  =  ^  -I-  -(b0  -ao)  H (b„  -  an),      i  =  1, ...,n-l. 

n  n 

Simple  Bounds 

In  many  applications,  it  is  helpful  to  introduce  simple  bounds  on  the  variables  [15]  such 
that  li  <  Xi  <  U{,  where  lj  and  Ui  are  lower  and  upper  bounds,  respectively.  This  is  a 
constrained  optimization  problem;  however,  it  has  a  particularly  simple  form.  A  solution 
x*  is  termed  feasible  if  it  satisfies  all  constraints.  A  step  5x  is  in  a  feasible  direction  if  it 
does  not  cause  x  to  exceed  any  of  the  bounds. 

The  following  discussion  treats  a  generic  iteration,  so  the  iteration  number  has  been  left 
out  for  simplicity.  At  any  iteration,  some  of  the  variables  will  be  equal  to  their  upper  or 
lower  bounds.  These  variables  xFX,  are  fixed  on  their  bounds  and  are  consequently  termed 
fixed  variables.  The  associated  bounds  are  placed  in  the  active  set.  The  remaining  free 
variables  xFR  are  allowed  to  vary  within  the  bounds.  The  search  direction  Sx  is  chosen  to 
change  only  the  free  variables.  Therefore,  6xFX  —  0  and  equation  (5.27)  reduces  to 

[HFR(x)]<5xFR  =  -gFR(x).  (5.28) 

With  5xFR  we  update  x  for  the  subsequent  iteration  as 


x  = 


xFR  +  dxFR 

XFX 
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If  5xFR  forces  xFR  to  violate  one  of  the  bounds,  that  bound  is  added  to  the  active  set  for 
the  subsequent  iteration. 

At  each  iteration,  the  active  set  must  be  analyzed  to  determine  if  the  constraints  in  the 
active  set  should  remain  active.  A  solution  x*  is  a  local  minimum  only  if  there  is  no  feasible 
step  8x  that  results  in  a  lower  value.  Thus,  if 

(<5x)Tg(x)  (5.29) 

is  non-zero  for  any  <5x,  then  x*  is  not  a  local  minimum.  Then  (5.29)  can  be  represented  as 

g(x).=  [A]TA 

where  [A]  is  the  set  of  constraint  coefficients  and  A  are  Lagrange  multipliers.  For  simple 
bounds,  [A]  is  composed  of  1  for  U  and  —1  for  u*,  so  the  Lagrange  multipliers  can  be 
estimated  as 

|        g(x)      lOTXi  =  k 

[  -gi(x)     for  Xi  =  ut 

Thus,  if  a  Lagrange  multiplier  for  an  active  simple  bound  is  negative,  the  corresponding 
bound  should  be  removed  from  the  active  set  for  the  subsequent  iteration. 


design  curve  [  design  curve 


(a)  Unconstrained  Cylindrical  Developable  Surface  (b)  Cylindrical  Developable  Surface 

With  Bounds 

Figure  5-3:  Effect  of  Simple  Bounds  on  Developable  Surfaces 

Many  solutions  are  possible  for  a  single  set  of  equations.  For  example,  Figure  5-3  shows 
two  equally  viable  solutions  to  the  same  set  of  equations.  Both  examples  have  design  curves 
with  control  points  (0,  0,  0),  (1.8,  3,  -0.7),  (3.3,  -2,  1.5)  and  (4,  0.5,  0)  and  endpoints  (0,  0, 
3)  and  (4,  0.5,  3)  for  the  second  directrix.  Figure  5-3(a)  is  unconstrained,  yielding  interior 
control  points  for  rg  of  (1.8,  3,  3)  and  (3.3,  -2,  3).  Figure  5-3(b)  is  unconstrained  in  the  x 
and  y  directions,  but  has  a  lower  bound  equal  to  5  set  in  the  z  direction.  This  produces  a 
resultant  curve  with  interior  control  points  (1.8,  3,  5)  and  (3.3,  -2,  5).  Since  both  cylinders 
are  developable  surfaces  with  C2  continuity  that  include  the  endpoints  specified  by  the  user, 
they  are  equally  valid  solutions  to  the  equations.  Both  surfaces  are  exactly  developable  with 
zero  Gaussian  curvature.  Although  this  example  is  rather  trivial,  it  demonstrates  the  control 
available  to  the  user  to  affect  the  shape  of  the  resulting  surface  by  specifying  bounds  on  the 
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solution.  Another  example  is  given  in  Section  5.3.2. 

Generalized  Elimination  Method 

From  equation  (5.23)  it  can  be  seen  that  Eq(x.)  =  0  reduces  to  a  linear  equation, 

(b0  -  ao)  x  (ai  -  ao)  •  (bi  -  b0)  =  0,  (5.31) 

where  bi  is  the  only  unknown.  This  condition  forces  the  first  four  control  points,  namely 
ao,  ai,  bo  and  bi,  to  lie  on  a  plane,  and  has  the  following  geometric  interpretation.  The 
two  Bezier  curves  r^u)  and  r^u)  pass  through  the  points  ao  and  bo  and  their  tangents 
are  in  the  directions  of  the  vectors  aoai  and  bobi.  Since  a  developable  surface  has  the  same 
tangent  plane  at  all  points  along  a  ruling  and  these  tangents  are  on  the  same  ruling,  u  =  0, 
they  must  be  coplanar.  Therefore,  ao,  ai,  bo  and  bi  must  lie  on  the  same  plane. 
Similarly,  the  last  equation,  i?3n_2(x)  =  0,  is  linear  in  bn_i: 

(b„  -  an)  x  (an  -  an_i)  •  (b„  -  b„-i)  =  0.  (5.32) 

The  linear  equations  (5.31)  and  (5.32)  can  be  used  to  reduce  the  dimensionality  of  the 
optimization  by  two.  The  new  objective  function  can  be  introduced  without  terms  Eq(x.) 
and  £fn_2(x),  i.e. 

F(x)  =  fij(x)  +  £22(x)  +  . . .  +  £32n-3 (*)• 

The  following  discusses  the  generalized  elimination  method  [13,  15],  in  which  the  number  of 
variables  in  the  optimization  is  reduced  by  the  number  of  independent  linear  constraints. 
In  general,  if  we  denote  the  number  of  variables  in  the  optimization  as  N  (i.e.  x  G  R^) 
and  the  number  of  linear  equations  as  M ,  the  set  of  linear  equations  can  be  written  as 

[A]x  =  b 

in  which  [A]  is  M  x  N  and  b  G  RM.  In  other  words,  the  zth  row  of  [A]  contains  the 
coefficients  corresponding  to  the  zth  linear  equation,  and  is  therefore  the  gradient  of  that 
equation. 

This  reduction  in  dimensionality  can  be  described  formally  in  terms  of  two  subspaces, 
which  can  be  described  as 

•  the  M-dimensional  subspace  defined  by  the  rows  of  [A]. 

•  the  N  —  M  complementary  subspace  of  vectors  orthogonal  to  the  rows  of  [A]. 

Let  [Y]  denote  any  N  x  M  matrix  whose  columns  form  a  basis  for  the  range  space  of  [A]T. 
In  fact,  [Y]  may  be  taken  as  [A]T.  Let  [Zj  denote  a.n  N  x  (N  —  M)  matrix  whose  columns 
form  a  basis  for  the  set  of  vectors  orthogonal  to  the  rows  of  [A],  so  that  [A][Z]  =  [0]. 
The  matrix  [Z]  has  linearly  independent  columns  zi,  Z2,...,  zn-m  which  are  in  the  null 
space  and  act  as  basis  vectors  for  the  null  space.  Since  [Y]  and  [Z]  define  complementary 
subspaces,  every  iV- vector  x  has  a  unique  expansion  as  a  linear  combination  of  the  columns 
of  [Y]  and  [Z] 

x  =  [Y]xy  +  [Z]x2,  (5.33) 
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where  M- vector  xy  is  the  range-space  portion  of  x  and  (N  -  M)-vector  x2  is  the  null-space 
portion  of  x.  By  using  the  chain  rule  on  the  derivative  of  equation  (5.33),  it  follows  that 

v..  =  [Z]TVI. 

Thus, 

VXzF(x  +  [Z]x)  =  [Z]TVxF(x)  =  [Z]Tg(x) 
is  the  reduced  gradient  vector  and 

V2XzF(x  +  [Z]x)  =  [Z}TV2xF(x)[Z]  =  [Z]T[H(x)][Z] 

is  the  reduced  Hessian  matrix.    Given  a  feasible  initial  point  xo  we  solve  a  reduced  size 

(N  -  M)  x  (N  -  M)  linear  system 

[Z]T[H(x)][Z]x2  =  -[Z]Tg(x) 
for  xz,  where  the  subsequent  iteration  is 

x  =  x+  [Z]xz. 


5.3.2     B-Spline  Developable  Surfaces 

In  this  section  we  merge  ra  Bezier  patches  together  by  joining  them  along  their  end  rulings 
with  C2  continuity  such  that  the  resulting  developable  surface  can  be  represented  by  a  single 
B-spline  surface.  The  design  procedure  is  the  same  as  that  of  the  Bezier  case,  except  that 
the  design  curve  is  a  B-spline  curve.  Since  multiplication  of  B-spline  functions  is  complex 
and  computationally  expensive,  we  first  extract  the  Bezier  segments  from  the  B-spline  curve 
by  knot  insertion  [5]  as  described  in  Section  3.3.  We  then  evaluate  equation  (5.23)  for  each 
Bezier  segment  and  add  the  resulting  equations  to  the  objective  function  after  squaring 
them.  If  equations  Eq  and  £m(3n-2)  are  used  m  tne  general  elimination  method,  they  are 
not  included  in  the  objective  function. 

Recall  from  Section  3.2  that,  for  a  degree  (n,  1)  B-spline  surface  with  m  Bezier  segments, 
the  C1  and  C2  continuity  conditions  at  the  joints  are  given  by 


hi(bni+i  -  bm)  =  hi+i(bni  -  bnj_i),       1  <  i  <  m  (5.34) 


and 


bnt-i  +  -7 —  (b„j_i  -  bnj_2)  =  bni+i  +  7-L-(bni+i  -  bni+2).  (5.35) 

hi  hi+i 

Equations  (5.34)  and  (5.35)  reduce  to  six  scalar  linear  equations  for  each  junction.  Using 
these  linear  equations  in  the  generalized  elimination  method  described  in  Section  5.3.1 
ensures  that  the  surfaces  will  be  exactly  C2  continuous;  otherwise  they  may  have  some 
residual  value  after  optimization  that  would  cause  discontinuity.  For  a  degree  (n,  1)  B-spline 
patch  with  m  Bezier  segments,  the  number  of  unknowns  and  linear  equations  is  3 (ran  —  1) 
and  6(ra  —  1)  -I-  2  respectively,  so  the  reduced  number  of  unknowns  is  3ra(n  —  2)  +  1. 

For  a  B-spline  surface,  the  initial  approximation  for  x  is  calculated  in  the  same  manner 
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as  for  the  Bezier  case,  expressed  as 

i  ran  —  i 

bt  =  aj  H (b0-a0)H (bmn  -  amn),      i  =  1,. . .  ,mn  -  1.  (5.36) 

ran  ran 

Once  the  control  points  of  the  second  directrix  are  computed,  any  knots  that  may  be 
removed  from  the  curve  while  leaving  the  curve  geometrically  and  parametrically  unchanged 
are  removed  using  the  knot  removal  technique  described  in  Section  3.3. 

Figure  5-4(a)  shows  the  design  curve  r^(if)  which  is  a  cubic  B-spline  curve  with  control 
points  (0,  4,  0),  (1,  2,  1),  (2,  4,  1),  (3,  1,  0)  and  (4,  2,  0)  and  knot  vector  (0,  0,  0,  0,  .5,  1, 
1,  1,  1).  Figure  5-4(b)  shows  the  design  curve  after  knot  insertion  with  the  resulting  two 
Bezier  segments  and  the  selected  endpoints  r#(0)  and  rs(l),  which  are  (0,  4,  3)  and  (8,  0,  3). 
When  this  problem  is  treated  as  an  unconstrained  optimization,  the  answer  converges  to  the 
trivial  solution  where  the  internal  control  points  of  the  new  curve  are  equal  to  the  internal 
control  points  of  the  design  curve  as  shown  in  Figure  5-4(c).  To  avoid  this  situation,  lower 
bounds  are  set  on  the  z  variables  where  li  —  3.0,  i  =  2,  5,  8,  11  and  14,  and  the  resulting 
internal  control  points  are  (2,  0,  5),  (4,  4,  5)  and  (6,  -2,  3).  The  resulting  surface  is  shown 
in  Figure  5-4(d)  with  two  Bezier  patches  and  the  control  mesh  delineated.  Figure  5-4(e) 
shows  the  final  developable  B-spline  surface  after  knot  removal. 

Aumann  [2]  proves  that  a  degree  (3-1)  Bezier  surface  whose  directrices  lie  in  parallel 
planes  is  developable  if  and  only  if  one  of  the  following  conditions  holds: 

1.  The  rulings  R(0,u)  and  R(l,v)  are  intersecting  or  parallel. 

2.  The  coefficient  of  the  cubic  term  of  the  design  curve  is  zero. 

3.  The  design  curve  contains  exactly  one  singular  point. 

In  constructing  examples  for  this  thesis,  it  was  found  that  skew  end  rulings  induce  non-zero 
Gaussian  curvature  into  (3-1)  Bezier  surfaces  whose  directrices  are  3D  space  curves  as  well. 

Algorithm 

The  following  steps  describe  the  algorithm  for  the  solution  method. 

1.  The  design  curve  r^(u)  and  the  endpoints  r^(0)  and  r#(l)  are  provided  by  the  user 
along  with  any  constraints  such  as  simple  bounds  or  constraining  planes  on  the  interior 
control  points  of  tb(u). 

2.  Bezier  segments  are  subdivided  out  of  the  B-spline  curve  by  knot  insertion  as  described 
in  Section  3.3. 

3.  Initial  values  for  x  =  (bg,  6g,  6q> . .  • ,  &™n,  ^mn^mn)  are  calculated  using  a  linear  inter- 
polation method  (5.36). 

4.  The  objective  function  is  optimized  in  two  steps.  First,  beginning  with  the  initial 
values,  an  initial  feasible  point  is  calculated  that  meets  all  constraints  and  bounds. 
Then,  an  optimal  solution  is  found  following  the  method  described  in  Section  5.3.1, 
using  the  following  steps: 

•  If  x  satisfies  the  objective  function  F(x)  =  0  and  the  Lagrangian  multiplier 
estimates  are  non-negative  for  all  active  bounds,  the  iteration  process  is  complete. 
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(a)  Design  Curve  r  A(u)  with  Control  Polygon  and  Knots 


(b)  Design  Curve  after  Knot  Insertion  with  Bezier  Segments 
and  Endpoints  of  r  g(u)  Shown 


(c)  Resulting  Surface  Without  Simple  Bounds 


(d)  Design  Curve  and  Resulting  Directrix  before  Knot  Removal 


(e)  Final  Developable  (3,1)  B-Spline  Surface 


Figure  5-4:  Design  of  Developable  B-Spline  Surfaces 
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(Note  that  the  objective  function  is  a  sum  of  squares,  so  the  minimum  possible 
sum  is  zero.) 

•  A  feasible  step  direction  <5x  is  found  using  (5.28). 

•  If  5x  causes  x  to  exceed  any  of  the  bounds,  those  bounds  are  added  to  the  active 
set. 

•  The  Lagrangian  multiplier  estimates  are  calculated  for  all  bounds  in  the  active 
set  using  (5.30).  If  any  are  negative,  the  associated  bound  is  removed  from  the 
active  set. 

•  x  is  set  to  x  +  <5x  and  the  next  iteration  is  begun. 

5.  Once  the  resultant  directrix  tb(u)  is  computed,  the  developable  B-spline  surface  is 
found  using  equation  (4.2). 

6.  Knots  that  can  be  removed  while  leaving  the  surface  unchanged  are  removed  using 
the  method  described  in  Section  3.2. 

7.  The  Gaussian  curvature  bound  is  calculated  using  the  method  described  in  Sec- 
tion 5.3.3  to  gauge  the  quality  of  the  resulting  surface. 

8.  The  surface  is  developed  onto  a  plane  using  the  technique  described  in  Section  5.4. 

5.3.3     Accuracy  of  Resulting  Developable  Surfaces 

Since  the  optimization  results  do  not  always  precisely  satisfy  F(x)  =  0,  we  need  a  measure 
to  determine  how  close  the  resulting  surface  is  to  an  exactly  developable  surface.  Recall 
from  Section  4.2  that  a  necessary  and  sufficient  condition  for  a  surface  to  be  developable  is 
for  the  Gaussian  curvature  to  be  zero.  This  property  can  be  used  as  a  measure.  For  each 
resulting  Bezier  patch  R(u,v)  =  YLf^oY^jzzO^UjBi^i^Bj^fa),  the  Gaussian  curvature,  K, 
can  be  expressed  in  terms  of  ratios  of  bivariate  Bernstein  polynomials  as  [31] 


K{u,v)     = 


(S  •  Ftuu)(S  •  Rw)  —  (S  ■  Huv) 

E^OAEUoSijBi,Bn-A(u)BM{v) 


(5.37) 


where  S  =  R^  x  R^  and  S  =  |Ru  x  R^.  If  R(w,u)  is  a  developable  surface,  then  KtJ  must 
be  zero  for  %  =  0, . . . ,  6n  —  4  and  j  =  0, . . . ,  2.  However,  for  surfaces  where  F(x)  /  0,  the 
Kij  coefficients  are  not  necessarily  always  precisely  zero. 
From  equation  (5.37),  we  have 


\K(u,v)\  = 


Knum{u,v) 


^■denomV^i  V) 


\Knum(u,v)\  (53g) 


\Kdenom(u,v)\ 


We  can  easily  find  the  upper  bound  of  \Knum(u,v)\  as 

6n-4    2 


\Knum{u,v)\     = 


^  ^2  KijBi>6n-4{u)Bjt2(v) 

i=0  j=0 
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6n-4    2 

6n-4    2 
<      WBx\Kij\    ^    ^2Bifin-4(u)Bj,2(v) 
13  i=0  j=0 

=     max|/fy|. 

The  denominator  Kdenom(u,v)  is  always  positive  if  the  surface  is  regular,  since  |Ru  x  R„| 
is  always  positive  for  a  regular  surface.  However,  some  of  the  coefficients  Sij  can  be  zero 
or  negative.  If  there  are  non-positive  Sij  we  can  subdivide  the  surface  into  smaller  patches 
until  all  Sij  become  positive.  This  is  always  possible  by  virtue  of  the  convergence  of  the 
control  polyhedron  to  the  surface  with  repeated  subdivision  [22].  Then  if  all  Sij  >  0  we 
have  the  inequality 


\Kdenom(u,v)\      = 


8n-4    4 

^2   ^2SijBifln-4{u)BjA{v) 
i=0  j-0 

8n-4    4 

>    min|5^'|  E  ^2Bi,8n-^(u)BiAv) 

13  i=0  j=0 

=     min  Sij . 
ij 

Therefore  we  can  rewrite  equation  (5.38)  as 

|JT(«..)I  <  ^^A.  (5.39) 

munj  bij 

If  minZJ  Sij  is  very  small,  the  bound  of  the  Gaussian  curvature  can  be  extremely  large. 
In  such  case  we  further  subdivide  the  patch  and  re-evaluate  the  bound  until  such  bound 
converges  to  a  stable  value. 

Alternatively,  the  curvature  may  be  analyzed  at  discrete  locations  on  the  surface,  which 
provides  a  rough  estimate  but  may  not  find  the  actual  bound  of  the  maximum  curvature. 

We  use  the  Gaussian  curvature  bound  to  give  a  physical  measure  of  the  developability  of 
the  surface.  Throughout  this  thesis  we  say  the  surface  is  exactly  developable  if  \K\  <  10-10. 
For  example,  the  curvature  of  the  surface  shown  in  Figure  5-4,  as  determined  using  ten  sub- 
patches  per  Bezier  patch  in  equation  (5.39),  is  zero,  so  the  surface  is  exactly  developable. 


5.3.4     Developable  Surfaces  Constrained  by  Planes 

In  many  engineering  applications,  one  or  both  of  the  directrices  are  restricted  to  planar 
curves  which  are  not  necessarily  parallel.  In  this  situation,  we  can  force  the  control  points 
b;  and  thus  the  directrix  r#(it)  to  lie  in  a  plane.  This  plane  is  determined  by  the  two  end 
points,  bo  and  bmn  for  a  degree  (n,  1)  B-spline  patch  with  m  Bezier  segments,  and  a  point 
g  chosen  by  the  user.  The  equation  of  the  plane  is  given  by 

ax  +  by  +  cz  +  d  =  0  (5.40) 
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where,  for  a  degree  (n,  1)  B-spline  patch  with  m  Bezier  segments, 


(%-9VWmn-9Z)-(bZ0-9ZWmn-9y) 

(bz0-9z)(bxmn-9x)-(bx0-9x)(bzmn-9z) 

{K-9xWmn-gy)-{bl-gywmn-gx) 

-agx  -bgy  -  cgz . 


For  each  control  point  bj  to  lie  on  the  plane,  its  x,  y  and  z  components  must  satisfy 
equation  (5.40).  These  linear  constraints  can  be  used  to  reduce  the  number  of  unknowns  in 
the  equation  using  the  generalized  elimination  method  described  in  Section  5.3.1. 

However,  if  the  (mn-1)  planar  constraints,  the  (6m— 6)  C1  and  C2  continuity  conditions 
and  the  two  linear  developability  conditions  E$  and  £'m(3n_2)  are  all  used  to  reduce  the 
number  of  variables,  the  problem  may  become  so  overconst rained  that  is  is  difficult  or 
impossible  to  find  a  feasible  solution.  For  example,  a  cubic  B-spline  design  curve  with 
two  Bezier  segments  has  15  unknowns  (bi,b2,b3,b4,b5)  and  13  linear  equations,  leaving 
just  2  variables  for  the  entire  problem.  In  this  case,  we  use  a  subset  of  the  total  available 
linear  equations  as  linear  constraints  and  add  the  remaining  equations  to  the  objective 
function  after  squaring  them.  In  selecting  the  linear  equations  to  be  used  as  constraints, 
it  is  important  to  remember  that  the  linear  constraints  may  be  satisfied  exactly  while 
the  objective  function  equations  may  not  be.  For  example,  if  C2  continuity  is  a  must, 
equations  (5.34)  and  (5.35)  are  used  as  linear  equality  constraints  and  the  equations  for 
planar  constraints  are  included  in  the  objective  function. 


(a)  Developable  Surface  with  Constraining  Plane 


(b)  Final  Degree  (3,1)  Developable  Surface 


Figure  5-5:  Planar  Constraint  on  Developable  B-Spline  Surface  Design 

Figure  5-5  shows  a  developable  B-spline  surface  that  was  designed  with  the  second 
directrix  constrained  to  a  plane.  The  same  design  curve  and  endpoints  were  used  as  those 
used  for  Figure  5-4,  except  the  second  directrix  was  constrained  to  lie  in  the  z  =  3  plane. 
Note  that  since  the  second  directrix  was  constrained  to  a  plane,  the  trivial  solution  was  not 
an  option,  and  simple  bounds  were  not  required.  The  continuity  conditions  were  used  as 
linear  constraints. 

5.3.5     Degenerate  Developable  Surfaces 

Sometimes  topologically  triangular  patches  are  required  to  produce  a  desired  surface  design. 
A  three  sided  patch  parameterized  over  a  rectangular  domain  can  be  constructed  by  allowing 
one  of  the  boundary  generators  to  be  zero  length,  say  R(0,i>).  Faux  and  Pratt  [11]  showed 
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that  the  surface  normal  at  this  point  may  be  evaluated  using  a  limiting  process  based  on  a 
Taylor  series  expansion 


lim  N(u,t>)  = 


u->0 


Ru(0,v)  xR„t(0,D) 
|Ru(0,u)  xRu^O.v)! 


(5.41) 


For  a  degree  (n-1)  developable  Bezier  patch,  equation  (5.41)  reduces  to 


N(0,v)  = 


(Pi,o-Po,o)x(Pi,i-Pi,q) 
l(Pi,o-Po,o)x(P1,1-P1,o)| 


Therefore,  the  unit  normal  vector  at  the  degenerate  point  always  exists  as  long  as  the  three 
control  points,  Po,o>  Pi,o  and  Pi,i,  do  not  coincide  with  each  other,  and  it  is  orthogonal  to 
the  plane  formed  by  the  three  control  points  as  shown  in  Figure  5-6. 

The  Gaussian  curvature  at  the  degenerate  point  cannot  be  calculated  using  the  method 
described  in  Section  5.3.3  because  the  denominator,  |Ru  x  R^l4,  vanishes.  Wolter  and 
Tuohy  [43]  devised  a  method  for  determining  Gaussian  curvature  at  a  degenerate  point  using 
the  second  order  partial  derivatives  of  a  local  height  function  representation  of  the  surface 
assuming  the  surface  has  a  tangent  plane  and  well  defined  curvatures  at  the  degenerate 
point.  For  this  description,  assume  that  the  degenerate  point  occurs  at  R(0,v).  The 
coordinate  system  is  shifted  such  that  the  origin  coincides  with  the  degenerate  point,  the 
z  axis  is  aligned  with  the  normal  to  the  surface  and  the  x  axis  is  aligned  with  the  tangent 
to  the  surface  in  the  non-degenerate  parameter  direction;  in  this  case,  Ru(0, 0).  The  y 
axis  is  orthogonal  to  the  others  to  form  a  right-handed  system.  Three  linearly  independent 
isoparameter  curves  in  u  are  selected  on  the  surface  such  that  an(u)  =  R(u, vn),  n  = 
1, . . . ,  3.  The  normal  curvature  in  the  direction 

0R(O,«n)       dan(0) 

™n  =  5 =  5 =  {Xn,yn) 


du 


du 


at  the  point  R(0, 0)  is  given  as 


n   kn    =    n-aJ(O) 


=    9xx(0, 0)a£  +  2gxy(0, 0)xnyn  +  gyyy\ 


which  can  be  represented  in  matrix  form  as 


where 


[Lla  =  b 


a     —      \9xxi9xyi9yy) 

b     =     (n-ai'(0),n-aJ(0),n-aS(0))r 


(5.42) 


and 


[Ll  = 


x\    2xi2/i    y\ 

A    2x22/2    v\ 
A     2x3y3    y\ 


If  the  fundamental  forms  are  converted  to  the  new  coordinate  system,  the  Gaussian  curva- 


65 


ture  can  be  represented  as 


-ft  —  9xx9yy       9xy 


(5.43) 


The  Gaussian  curvature  at  a  degenerate  point  can  thus  be  determined  by  solving  (5.42)  for 
a  and  substituting  the  resulting  values  for  gxx,  gxy  and  gyy  into  (5.43). 

To  estimate  the  maximum  absolute  value  of  the  Gaussian  curvature  bound  of  a  degen- 
erate surface  with  the  degenerate  point  at  R(0,0),  the  surface  is  first  subdivided  at  u  =  e, 
where  e  is  a  small  positive  number  so  that  the  resulting  patch  with  e  <  u  <  1  is  a  nonde- 
generate  patch.  The  Gaussian  curvature  is  calculated  at  u  =  0  using  equation  (5.43),  and 
the  b.ound  for  the  Gaussian  curvature  for  the  nondegenerate  patch  is  evaluated  using  the 
technique  described  in  Section  5.3.3. 
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Figure  5-6:  Triangular  Degenerate  Cubic  Bezier  Patch 


Figure  5-6  shows  a  degenerate  cubic-linear  Bezier  patch  with  non-planar  directrices.  The 
design  curve  is  made  up  of  the  control  points  (0,  0,  0),  (1.8,  3,  -0.7),  (3.3,  -2,  1.5)  and  (4, 
0.5,  0).  The  selected  end  points  are  (0,  0,  0)  and  (4,  0.5,  2).  The  Gaussian  curvature  at  the 
degenerate  point  was  found  to  be  8.91  x  10" 14.  The  surface  was  split  at  the  parameter  value 
u  =  0.0001  resulting  in  a  maximum  bound  for  the  Gaussian  curvature  of  zero.  Therefore, 
the  surface  is  exactly  developable.  Note  that  there  are  only  three  points  which  define  the 
end  rulings  and  they  always  lie  on  the  same  plane.  Therefore,  there  is  no  skew  between  the 
two  end  rulings  and  the  surface  often  results  in  an  exact  developable  surface. 


5.4     Development  of  a  Developable  Surface  onto  a  Plane 

In  the  manufacture  of  developable  surfaces,  it  is  necessary  to  calculate  the  plane  development 
of  these  surfaces.  This  determines  cut  information  for  the  flat  piece  of  metal  so  it  can  be 
cut  to  the  proper  dimensions  before  it  is  rolled. 

This  formulation  follows  the  method  of  Faux  and  Pratt  [11,  p. 274],  using  the  fact  that  a 
developable  surface  may  be  mapped  isometrically  onto  a  plane.  The  Frenet-Serret  formulae 
state  that 


t'  =  Kn 


(5.44) 
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where  t  is  the  unit  tangent  vector  to  a  curve,  n  is  the  unit  normal  vector  to  a  curve  and 
k  is  the  curvature.  For  a  planar  curve  in  the  (X,  Y)  plane,  we  can  define  the  unit  normal 
vector  as 

n  =  ez  x  t  (5.45) 

where  ez  =  (0,0, 1).  Substituting  (5.45)  into  (5.44)  yields 

d2X        dY      n 

—j-  -  k—  =  0  (5.46) 

as*  as 

d2Y        dX      n 

+  K        =0  (5.47) 

asz  as 

where  (X,  Y)  denote  the  2D  coordinates  on  the  developed  plane  (X,Y).  If  we  rewrite  equa- 
tion (5.46)  in  terms  of  the  parameter  u,  we  obtain 

d2X        dY  d  (dX\         dY 

—  K—r-      =      —      -T—      —  K 


ds2  ds  ds  \  ds  J         ds 

d   t  dX  du\         dY  du 
ds  \du  ds  J         du  ds 


d2X  fdu\2      dX  (dPu\         d^du_ 


du2    \ds  J         du   \  ds2  I  du  ds 


If  (du/ds)2  7^  0,  this  reduces  to 


d2X       (0)  dX 


k     dY 


du2        ( du\2  du        (du\  du 
\ds)  \dsJ 

Similarly,  equation  (5.47)  can  be  represented  as 


+  T^-T--7TX-T-=0-  (5-48) 


A^M"    •     «    «.a  (5.49) 


du2  (  du\2  du  (  du\ 

\ds)  V  ds  ) 


du 


The  development  is  based  on  the  fact  that  curves  on  isometric  surfaces  have  the  same 
geodesic  curvatures.  It  can  be  seen  from  Figure  2-4  that  for  a  planar  curve,  since  the  angle 
between  the  surface  normal  and  the  normal  to  a  curve  is  90  degrees,  the  curvature  vector 
equals  the  geodesic  curvature  vector.  Therefore,  k  in  equations  (5.48)  and  (5.49)  can  be 
replaced  by  k9,  the  geodesic  curvature  of  the  curve  on  the  developable  surface.  If  we  choose 
the  curve  on  the  developable  surface  to  be  an  isoparametric  curve  in  terms  of  u,  R(u,vn), 
we  can  replace  du/ds  and  d2u/ds2  in  equations  (5.48)  and  (5.49)  by 

du  1 

—     =     :=— ; tt         and 

ds  |Ru(u,Vn)| 

d2u  Ru{u,vn)  ■  Kuuju.Vn) 

ds2  Ru(w,un)  •  Ru(u,vn)  ' 
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Thus  we  have 


dX 

du 

=    P 

dY 

du 

=     Q 

dp 

(.-Tt-tl    '   *^Uu) 

1  "D 

du 

(R-u  -  R-u) 

—  qKgitVu 

dq 

^■*MI    '   *^UU  J 

1  T"» 

du 

(Ru  •  ru) 

+  PKg\tlu 

(5.50) 

The  surface  will  be  developed  onto  the  plane  oriented  such  that  the  point  R(0,0)  is 
located  at  the  origin  of  the  (X,  Y)  coordinate  system  in  the  plane  and  the  vector  Ru(0,0) 
coincides  with  the  X  axis,  as  shown  in  Figure  5-7.  We  integrate  the  system  (5.50)  along 
the  directrix  that  corresponds  to  v  =  0  from  u  =  0  to  u  =  1,  using,  for  example,  a  Runge 
Kutta  scheme.  Since  isometric  maps  are  conformal,  the  angle  between  the  directrix  and  the 
generator  at  (0, 0)  is  the  same  in  both  representations  and  can  be  found  by 

RuCCO)       1^(0,0) 
cost/  = 


|Ru(0,0)|     |Rt,(0,0)| 

Straight  lines  on  the  surface  map  to  straight  lines  on  the  plane,  so  the  generator  length 
on  the  surface  is  measured  and  transferred  to  the  plane,  giving  the  initial  point  for  the  de- 
velopment of  the  second  directrix.  Then  the  system  of  ordinary  differential  equations  (5.50) 
is  integrated  using  a  fourth  order  Runge-Kutta  method  [38]  along  the  directrix  that  cor- 
responds to  v  =  1  from  u  =  0  to  u  =  1.  Finally,  the  endpoints  of  the  two  directrices  are 
connected. 

To  verify  that  the  development  of  the  surface  is  accurate,  one  can  check  the  results 
by  computing  the  angle  between  the  directrix  corresponding  to  v  =  0  and  the  generator 
corresponding  to  u  =  1,  computing  the  length  of  that  generator,  and  determining  the 
location  of  the  point  R(l,l)  in  the  plane.  Comparing  this  point  to  the  point  calculated 
above  determines  the  accuracy  of  the  development.  The  development  will  be  more  accurate 
with  increased  steps  in  the  Runge-Kutta  scheme. 

The  surface  that  was  designed  in  Section  5.3.2  is  developed  onto  a  plane  here.  This 
surface  is  shown  in  Figure  5-7(a).  It  was  developed  onto  a  plane  using  31  steps  in  a  fourth 
order  Runge-Kutta  routine  as  shown  in  Figure  5-7(b).  The  accuracy  for  the  check  at  R(l,  1) 
was  10-6. 

5.5     Examples 

5.5.1     Example  of  a  Developable  B-Spline  Strip  Surface 

This  example  is  a  windshield  of  an  automobile  as  depicted  in  Figure  5-8.  Since  automobile 
windshields  are  generally  sandwiched  between  two  parallel  planes,  we  can  use  the  technique 
introduced  in  Section  5.2.  The  design  curve,  r^(ti),  is  given  in  the  xy-plane  by  a  cubic 
B-spline  curve  with  control  points  (0,  2.1,  0),  (0.5,  0.9,  0),  (1.5,  0.3,  0),  (4.5,  0,  0),  (7.5,  0.3, 
0),  (8.5,  0.9,  0),  (9,  2.1,  0)  and  a  knot  vector  (0,  0,  0,  0,  0.3,  0.5,  0.7,  1,  1,  1,  1).  The  two 
end  boundary  conditions  are  chosen  to  be  bo  =  (0.2,  3.5,  1.5)  and  b6  =  (8.8,  3.5,  1.5).  The 
design  curve  is  first  split  into  four  cubic  Bezier  curves  and  the  control  points  of  r#(it)  are 

68 


R(O.O) 


directrix,  rB(u) 


directrix,  rA(u) 


(a)  Developable  B-Spline  Surface 


(b)  Plane  Development  of  Surface 


Figure  5-7:  Development  of  B-Spline  Developable  Surface 

evaluated.  The  resulting  B-spline  curve,  rjg(u),  has  internal  control  points  (0.56,  2.64,  1.5), 
(0.96,  2.18,  1.5),  (1.84,  1.81,  1.5),  (3.27,  1.62,  1.5),  (4.5,  1.53,  1.5),  (5.73,  1.62,  1.5),  (7.16, 
1.81,  1.5),  (8.04,  2.18.  1.5)  and  (8.44,  2.64,  1.5),  with  knot  vector  (0  0  0  0  0  .3  .3  .5  .5  .7 
.711111).  Note  that  the  resulting  surface  is  not  cylindrical.  After  degree  elevating  the 
design  curve  by  one  and  removing  the  knots  from  r#(tt),  we  obtain  the  degree  (4-1)  B-spline 
developable  surface  shown  in  Figure  5-8.  The  windshield  surface  is  then  developed  onto  the 
plane  as  shown  in  Figure  5-8(c). 

5.5.2     Example  of  a  Developable  B-Spline  Surface  with  3D  Directrices 

The  stack  of  a  ship,  depicted  in  Figure  5-9,  was  constructed  using  a  B-spline  design  curve 
with  the  resulting  curve  constrained  to  a  plane  not  parallel  to  the  plane  of  the  design  curve. 
Therefore,  the  method  described  in  Section  5.3  was  required.  The  stack  was  split  along  the 
plane  of  symmetry  (the  xz  plane)  from  forward  to  aft.  The  simple  shape  allowed  definition 
of  the  design  curve  with  a  single  B-spline  curve  of  two  Bezier  segments  for  each  half  stack. 
A  cubic  design  curve  was  selected  that  defined  the  base  of  the  stack,  with  control  points  (0, 
0,  0),  (0,  10,  0),  (15,  10,  0),  (35,  3,  0)  and  (35,  0,  0)  and  knot  vector  (0,  0,  0,  0,  .5,  1,  1,  1, 
1).  The  two  endpoints  were  selected  to  be  bo  =  (6,  0,  31)  and  b6  =  (24,  0,  29).  The  second 
directrix  was  constrained  to  lie  on  a  plane,  with  the  third  point  to  define  the  plane  being 
(24,  10,  29).  The  optimization  process  was  conducted  using  equations  (5.34)  and  (5.35)  to 
reduce  the  number  of  variables,  and  no  simple  bounds  were  employed.  The  resulting  surface 
was  degree  (3,1).  The  maximum  Gaussian  curvature  determined  using  equation  (5.39)  with 
10  Bezier  subpatches  was  2.19  x  10~7.  The  surface  was  then  reflected  with  respect  to  the 
xz  plane,  resulting  in  a  symmetric  body.  One  of  the  surfaces  surfaces  was  developed  onto 
a  plane,  as  shown  in  Figure  5-9  (c). 
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(b)  Final  Developable  Surface 


(c)  Surface  Developed  onto  Plane 


Figure  5-8:  Automobile  Windshield,  a  Developable  Degree  (4-1)  B-spline  Strip  Surface 
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(a)  Plan  View  with  Control  Points 


(b)  Sheer  View 


(c)  Half  Stack  Developed  onto  Plane 


design  curve 


Figure  5-9:  Ship  Stack,  A  Developable  Degree  (3,1)  B-Spline  Surface  with  \K\  <  2.19  x  10"7 
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Chapter  6 

Geodesies  on  Developable  Surfaces 


6.1  Introduction 

Geodesies  axe  useful  in  ship  design  and  shipbuilding  for  determining  the  layout  of  seams  and 
butts  in  a  ship  hull,  as  described  by  Munchmeyer  and  Haw  [33].  This  chapter  shows  that  a 
two  point  boundary  value  problem  (BVP)  for  computing  geodesies  on  a  developable  surface 
can  always  be  reduced  to  an  initial  value  problem  (IVP),  with  a  corresponding  savings  in 
time  for  the  solution  process.  Most  problems  that  arise  in  applications  of  geodesies  are 
BVPs,  which  are  much  more  difficult  to  solve  than  IVPs. 

6.2  Formulation 

Recall  from  Section  2.3  that  the  governing  equations  of  a  geodesic  are  given  by  a  set  of 
coupled  second  order  ordinary  differential  equations  (2.29)  and  (2.30) 


d2u      „,    / du\2      „,„!   dudv      _i    ( dv\2 

d2v  ,r2    (<*u\2         2  dudv        2    (dv\2 
d?  +  in\Ts)    +2il2d~sd~s+l22{d~s)    -° 


where  u  and  v  are  related  by  the  first  fundamental  form  ds2  =  Edu2  +  2Fdudv  +  Gdv2. 
These  two  second  order  differential  equations  can  be  rewritten  as  a  system  of  four  first  order 
differential  equations 


du 

ds 
dv 

ds 

4   =   -r!lP2-2r!2Pg-r^2  (6.3) 

dq 
ds 


=    P  (6.1) 

=    q  (6.2) 

!lP2-2r}2Pg-r22g2 

^        -r?lP2  -  2r22P<7  -  r2292.  (6.4) 


We  can  solve  the  system  as  an  initial- value  problem  (IVP),  where  all  four  boundary 
conditions  are  given  at  one  point,  or  as  a  boundary- value  problem  (BVP),  where  the  four 
boundary  conditions  are  specified  at  two  distinct  points.   Most  of  the  problems  that  arise 
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in  applications  of  geodesies  are  BVPs,  which  are  much  more  difficult  to  solve  than  IVPs. 

The  solution  of  an  IVP  is  unique;  however  for  a  BVP  it  is  possible  that  the  differential 
equations  will  have  many  solutions  or  even  no  solution  [25].  General  methods  for  the 
solutions  of  two-point  BVPs  can  be  found  in  [12,  25].  There  are  two  commonly  used 
approaches  to  the  numerical  solution  of  a  BVP,  namely  the  shooting  method  and  the  finite 
difference  method  (relaxation  method).  The  shooting  method  was  originally  used  to  adjust 
the  settings  of  artillery  equipment  to  hit  a  target.  To  implement  the  method,  the  unknown 
boundary  conditions  are  guessed  at  the  initial  point  and  an  IVP  method  is  implemented. 
Based  on  the  resulting  values  at  the  end  point,  the  initial  guess  is  adjusted  until  the  "target" 
is  "hit".  The  second  approach  is  based  on  a  finite  difference  approximation  to  -^  where 
y  =  (u,v,p,  q)  on  a  mesh  of  points  in  the  interval  [A,  B].  This  method  starts  with  an 
initial  guess  and  improves  the  solution  iteratively.  Maekawa  [30]  solves  the  geodesic  BVP 
on  free-form  parametric  surfaces  by  the  finite  difference  method.  In  general,  the  shooting 
method  is  very  sensitive  to  the  initial  guess  at  point  A,  and  is  much  less  stable  than  the 
finite  difference  method  [30]. 

For  a  B-spline  developable  surface  patch  there  is  a  single  unique  solution  to  the  system. 
Since  a  geodesic  on  a  developable  maps  to  a  straight  line  on  a  developed  plane,  there  is 
only  one  solution  to  the  system  on  a  developable  surface.  Here  we  exclude  periodic  surfaces 
such  as  cylinders  where  there  can  be  more  than  one  solution,  as  shown  in  Figure  6-1. 


Figure  6-1:  Geodesies  on  a  Closed  Surface 

In  this  section  it  is  shown  that  all  two  point  BVPs  for  solving  geodesies  on  developables 
can  be  reduced  to  IVPs  using  the  following  properties  described  in  Section  4.2: 

•  A  developable  surface  can  be  mapped  isometrically  onto  a  plane. 

•  Isometric  surfaces  have  the  same  Gaussian  curvature  at  corresponding  points. 

•  Corresponding  curves  on  these  surfaces  have  the  same  geodesic  curvature  at  corre- 
sponding points. 

•  Every  isometric  mapping  is  conformal;  i.e.  the  angle  of  intersection  of  every  arbi- 
trary pair  of  intersecting  arcs  on  a  developable  surface  is  the  same  as  that  of  the 
corresponding  inverse  image  in  the  plane  at  the  corresponding  points. 
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•  A  geodesic  on  a  developable  surface  maps  to  a  straight  line  in  the  plane. 

The  basic  procedure  is  to  map  the  two  desired  points  on  the  developable  surface  to  a  plane, 
draw  the  straight  line  between  them  and  determine  the  angle  between  the  generator  and 
the  geodesic  line  at  one  of  the  end  points.  The  angle  can  be  used  to  determine  du/ds  and 
dv/ds.  Thus,  all  the  information  required  for  an  IVP  is  available. 


(a)  Surface  Developed  onto  Plane 


(b)  uv  Parametric  Space 


R(0,0) 


(c)  Geodesic  on  3D  Surface 


(d)  Final  3D  Surface  with  Geodesic 


Figure  6-2:  Geodesic  on  a  Degree  (3,1)  Developable  Surface 

Given  two  points  A  and  B  on  the  developable  surface  R(u,  v)  as  shown  in  Figure  6-2(c), 
the  corresponding  points  {Xa,  Ya)  and  (Xb,  Yb)  in  the  developed  plane  are  required.  To  find 
the  points  in  the  plane,  (Xq,  Yq)  is  set  as  the  (0, 0)  point  in  the  plane  corresponding  to  R(0, 0) 
on  the  surface.  Using  the  method  presented  in  Section  5.4,  the  directrix  corresponding  to 
v  =  0  is  developed  into  the  plane  to  determine  the  point  C  =  (ua,0),  shown  in  Figure  6- 
2(a).  The  angle  between  the  directrix  and  the  generator  at  (ua,0)  is  the  same  in  both 
representations  and  can  be  found  by 


cos#  = 


Ru{uA,0)       Cl 
|Ru(tM,0)|     i^Ji 


where  CA  is  an  isoparametric  line  R(u^,  v)  which  is  a  straight  line  on  the  surface.  Therefore 
it  is  a  geodesic  and  will  be  developed  into  the  plane  as  a  straight  line.  The  distance  \CA\ 
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is  given  by 

\Cl\  =  \J{XA  ~  XC)2  +  (VA  ~  VC)2  +  {ZA  ~  ZC)2. 

The  point  A  on  the  plane  is  found  using  C,  \CA\  and  9.  The  point  B  on  the  plane  is  found 
by  following  the  same  procedure,  and  the  points  are  connected  as  shown  in  Figure  6-2 (a). 

The  angle  u  between  CA  and  AB  is  given  by 

AB*-c2 
cos  a;  =  — . 

\AB\\CA\ 

This  angle  u>  is  preserved  between  the  isoparametric  line  t(ua,v)  and  the  geodesic  curve 
g(s)  on  the  developable  surface  at  point  A.  Thus  we  have 

cosu;  = 


IR.IIg'WI 

where  the  tangent  vector  to  the  geodesic  is  given  by 

..  .      _    du      _   dv  ,„  „x 

g'(s)  =  Ku—  +  Kv  —  .  (6.5) 

as  as 

The  angles  u  and  6  are  shown  in  Figures  6-2(a)  and  (c).  Multiplying  (6.5)  by  Ry  yields 

du  dv 

R,  •  g'(s)  =  R„   R„  —  +R„  •  R,,—  =  cos(o;)|R„||g'(5)| 

ds  ds 

which  (since  |g'(s)|  =  1)  can  be  reduced  to 

„du      „dv  .   .    /^= 

F—  +  G—  =  cos(co)VG 
ds  ds 

where  F  =  Ru  •  R^  and  G  =  Ky    R„  (coefficients  of  the  first  fundamental  form).  Thus, 

dv       cos(cj)       F  (du\ 

(6.6) 


From  the  first  fundamental  form, 

//  »      //  v       r,  f  du\        ^^dudv       „  { dv\  .„  ... 

*(.)•*(.) -*(s)    +2F--+G(-)    =1.  (6.7) 

Plugging  (6.6)  into  (6.7)  and  solving  for  du/ds  yields 


du  _        js'm2(uj)G 


ds  V  EG-  F2 

and  thus  (6.6)  reduces  to 

dv       cos(u;)    ,1  F  .    .   . 

—  =  — t=^  ±  — =    . =  sin(a>). 

ds         y/G        VG  VEG  -  F2 
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Hence,  we  have  all  the  initial  conditions  required  to  solve  the  IVP  (equations  (6.1)  to  (6.4)) 
for  a  geodesic.  The  solution  to  the  IVP  yields  the  uv  parametric  values  for  the  geodesic  that 
are  graphed  in  Figure  6-2(b).  The  corresponding  three-dimensional  coordinate  values  are 
shown  in  Figures  6-2(c)  and  (d).  The  developable  degree  (3,1)  surface  shown  in  Figure  6-2 
has  the  control  points  and  knot  vectors  shown  in  Table  6.1.  The  geodesic  runs  from 

(uA,vA)  =  (0.1, 0.3)     to     (uB,  vB)  =  (0.9, 0.8). 


Po,o 

(-2,1,0) 

Po.i 

(  0,  2,  -3) 

Pi,o 

(0,-1,-1) 

Pi,o 

(1,  1,  -3.5) 

P2,0 

(2,  5,-1.4) 

P2,o 

(2,  4,  -3.7) 

P.3,0 

(4,  1,  -1) 

P3,0 

(3,  2,  -3.5) 

P4.0 

(6,  5,  0) 

P4,0 

(4,  4,  -3) 

Tu 

0,  0,  0,  0, 
•5,  1,1,1,1 

Tv 

0,  0,  1,  1 

Table  6.1:  Control  Points  for  a  Degree  (3,1)  Developable  Surface 


6.3     Examples 

In  this  section  we  conduct  numerical  experiments  to  compare  the  performance  of  the  IVP  to 
the  finite  difference  based  BVP  on  developable  surfaces.  For  a  fair  comparison,  we  determine 
the  minimum  number  of  steps  in  the  BVP  to  reliably  meet  the  given  tolerance  e  =  10-5 
for  Newton's  method,  see  [31].  We  then  find  the  minimum  number  of  steps  for  the  IVP 
to  achieve  the  same  order  of  accuracy  at  the  endpoint  B.  The  examples  were  run  on  a  180 
MHz  Silicon  Graphics  workstation. 

The  example  shown  in  Figure  6-2  runs  using  an  IVP  solution  with  66  steps  in  0.25 
seconds,  as  compared  to  0.60  seconds  for  a  BVP  solution  with  135  steps  and  7  iterations. 
This  is  over  twice  as  fast,  with  a  very  simple  surface  and  a  geodesic  that  is  very  close  to  a 
straight  line  in  uv  space. 

A  more  complex  degree  (4,1)  developable  surface  with  five  Bezier  patches  also  runs  over 
twice  as  fast:  0.88  seconds  for  the  IVP  solution  with  194  steps  instead  of  2.28  seconds  for 
the  BVP  solution  with  294  steps  and  11  iterations.  The  control  points  and  knot  vectors  for 
this  example  are  displayed  in  Table  6.2.  The  geodesic  runs  from  (ua,va)  =  (0.1,  0.9)  to 
{ub,vb)  =  (0.9,  0.1).  The  surface  and  its  geodesic  are  displayed  in  Figure  6-3. 

The  final  example  is  a  handkerchief- like  surface.  This  degree  (4,1)  developable  surface 
with  four  Bezier  patches  runs  in  0.50  seconds  for  the  IVP  and  2.66  seconds  for  the  BVP, 
which  is  over  5  times  faster.  The  IVP  required  176  steps,  while  the  BVP  required  320 
steps  and  12  iterations.  The  control  points  and  knot  vectors  for  this  example  are  shown  in 
Table  6.3.  The  geodesic  runs  from  (ua,va)  =  (0,0)  to  {ub,vb)  =  (1,1)-  The  surface  and 
its  geodesic  are  displayed  in  Figure  6-4. 
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(a)  Surface  Developed  onto  Plane 


(b)  uv  Parametric  Space 


(c)  Geodesic  on  3D  Surface 
Figure  6-3:  Geodesic  on  a  Degree  (4,1)  Developable  Surface 


Po,o 

(0.00,  3.00,  0.00) 

Po,i 

(0.00,  3.00,  -5.00) 

Pi.o 

(0.75,  1.50,  0.00) 

Pi,i 

(2.42,-1.84,  -5.00) 

P2,0 

(1.25,  1.50,  0.00) 

P2,l 

(3.88,-1.53,  -5.00) 

P3,0 

(1.96,  2.58,  0.00) 

Pa,i 

(5.85,  1.34,  -5.00) 

P4,0 

(2.50,  2.00,  0.00) 

P4,l 

(7.12,-2.28,  -5.00) 

P5,o 

(3.00,  1.17,  0.00) 

P5,l 

(8.21,-1.97,  -5.00) 

P6,0 

(3.50,  1.00,  0.00) 

Pe,i 

(9.09,-2.20,  -5.00) 

P7,0 

(4.00,  1.08,  0.00) 

Pr,i 

(9.86,-2.03,  -5.00) 

P8,0 

(4.50,  1.50,  0.00) 

Ps,i 

(1.04,-1.53,  -5.00) 

P9,o 

(5.04,  2.29,  0.00) 

P9,l 

(1.09,-6.81,  -5.00) 

Pio,o 

(5.75,  4.25,  0.00) 

PlO,l 

(1.13,  4.04,  -5.00) 

Pn,o 

(6.75,  4.50,  0.00) 

Pll,l 

(1.18,  1.05,  -5.00) 

Pl2,0 

(9.00,  3.00,  0.00) 

Pl2,l 

(1.20,-1.11,  -5.00) 

Tu 

0,  0,  0,  0,  0,  .2,  .2, 

.4,  .4,  .6,  .6,  .8,  .8, 

1,  1,  1,  1,  1 

Tv 

0,  0,  1,  1 

Table  6.2:  Control  Points  for  a  Degree  (4,1)  Developable  Surface 
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(a)  Surface  Developed  onto  Plane 


(b)  uv  Parametric  Space 


(c)  Geodesic  on  3D  Surface 
Figure  6-4:  Geodesic  on  a  Degree  (4,1)  Handkerchief- like  Developable  Surface 


Po,o 

(0.50,  3.00,  0.00) 

Po.i 

(0.50,  3.00,  -5.00) 

Pi,o 

(0.88,  1.50,  0.00) 

Pi,o 

(0.88,  1.50,  -5.00) 

P2,0 

(1.25,  2.00,  0.00) 

P2,0 

(1.25,  2.00,  -5.00) 

P3,0 

(1.96,  4.17,  0.00) 

P3,o 

(1.96,  4.17,  -5.00) 

P-4,0 

(2.50,  3.00,  0.00) 

P4,0 

(2.50,  3.00,  -5.00) 

P5,o 

(3.00,  1.67,  0.00) 

P5,o 

(3.00,  1.67,  -5.00) 

Pe,o 

(3.50,  3.00,  0.00) 

Pe,o 

(3.50,  3.00,  -5.00) 

Pt,o 

(4.04,  4.17,  0.00) 

Pt,o 

(4.04,  4.17,  -5.00) 

P8,o 

(4.75,  2.00,  0.00) 

Ps,o 

(4.75,  2.00,  -5.00) 

P9,0 

(5.13,  1.50,  0.00) 

P9,0 

(5.13,  1.50,  -5.00) 

Pio,o 

(5.50,  3.00,  0.00) 

Pio,o 

(5.50,  3.00,  -5.00) 

Tu 

0,  0,  0,  0,  0, 

.25,  .25,  .5,  .5,  .75,  .75, 

1,  1,  1,  1,  1 

Tv 

0,  0,  1,  1 

Table  6.3:  Control  Points  for  a  Degree  (4,1)  Handkerchief-like  Developable  Surface 
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Chapter  7 

Lines  of  Curvature  on  Developable 
Surfaces 

7.1  Introduction 

This  chapter  discusses  the  behavior  of  the  lines  of  curvature  on  a  developable  surface  near 
the  inflection  line.  A  developable  surface  does  not  possess  isolated  generic  flat  points  but 
rather  may  contain  a  line  of  non-generic  flat  points  along  the  generator.  This  property  can 
be  used  to  detect  an  inflection  line  where  the  (only)  non-zero  principal  curvature  whose 
principal  direction  is  perpendicular  to  the  generator  changes  its  sign.  Munchmeyer  [33] 
indicates  that  a  developable  surface  can  be  shaped  purely  by  rolling  and  should  be  fed  into 
the  roller  so  that  the  direction  of  the  generator  is  parallel  to  the  rolls.  However,  when  the 
non-zero  principal  curvature  changes  sign,  the  plate  can  no  longer  be  fed  into  the  roller 
in  the  same  direction,  since  the  bending  direction  changes  at  the  inflection  line.  Locating 
the  inflection  line  prior  to  fabrication  allows  the  manufacturer  to  either  cut  the  plate  along 
the  inflection  line  or  reverse  the  direction  the  plate  is  fed  through  the  rollers  at  that  line. 
Surface  inflection  was  studied  by  Hoitsma  [21],  who  showed  that  a  surface  has  an  inflection 
at  a  point  P  if  and  only  if  its  mean  curvature  changes  sign  in  the  neighborhood  of  P.  This 
chapter  further  extends  this  result. 

7.2  Inflection  Line 

This  section  establishes  the  definition  of  an  inflection  line  on  a  developable  surface.  First, 
some  additional  differential  geometry  properties  are  required. 

Recall  that  throughout  this  thesis,  it  is  assumed  that  the  u  =  const  isoparametric  line 
corresponds  to  the  generator  of  the  developable  surface  or,  in  other  words,  the  straight  line 
ruling  is  in  the  v  direction.  With  this  assumption,  R^  =  0;  consequently,  from  (2.13), 
N  =  R^  •  N  =  0.  From  (2.23),  since  Gaussian  curvature  of  a  developable  surface  is  zero, 

LN-M*_      -M2 
K~  EG-F*  -EG^F*-0  {7A) 

and  we  have  M  =  0.  Therefore,  mean  curvature  as  expressed  by  equation  (2.24)  reduces  to 

2FM-EN-GL  -GL 

2{EG-F2)  2{EG-F2)'  [      ' 
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Recall  from  (4.22)  that  the  nonzero  principal  curvature  is  given  by  k*  =  2H.  When  L  =  0, 
H  and  k*  become  zero;  otherwise,  k*  =  2H  ^  0. 

Next,  we  will  show  that  the  u  =  const  parametric  straight  lines  become  the  lines  of  zero 
curvature.  This  can  be  seen  from  the  fact  that 

•  the  u  =  const  parametric  straight  lines  have  zero  normal  curvature,  and 

•  no  other  direction  has  zero  normal  curvature. 
The  second  fact  comes  from  Euler's  theorem  (2.25), 

nn  =  «i  cos  (a)  +  «2  sin  (a). 

When  «2  =  0,  «n  =  «i  cos2 (a),  which  becomes  zero  only  when  a  =  \  or  ^,  corresponding 
to  the  direction  of  «2-  Similarly,  when  K\  =  0,  nn  =  0  only  when  a  =  0  or  7r,  corresponding 
to  the  direction  of  K\. 

Theorem  7.2.1  A  developable  surface  does  not  possess  generic  isolated  flat  points  1  but 
rather  may  contain  a  line  of  non-generic  flat  points  along  a  generator. 

Proof:  From  equation  (7.2),  L  vanishes  at  a  flat  point  (u/,Vf)  where  K  —  H  =  0.  There- 
fore, from  equation  (2.14)  we  have 

L{uf, vf)  =  -Hu(uf, Vf)  ■  Nu{uf, vf)  =  0.  (7.3) 

From  equation  (7.1),  M  =  0  on  a  developable  surface.  Therefore,  from  (2.15)  we  have 

M (uf, vf)  =  -R„(ti/, vf)  ■  Nu(u/, vf)  =  0.  (7.4) 

Since  N  is  a  unit  vector,  we  also  have 

N(iif,vf)-Nu(uf,vf)=0.  (7.5) 

If  Nu  is  not  zero,  then  from  (7.3),  (7.4)  and  (7.5)  Nu  must  be  perpendicular  to  R^,  R^ 
and  N.  This  is  impossible  because  N  is  perpendicular  to  both  Ru  and  R^,  and  R^  is  not 
parallel  to  R^.  Thus,  Nu(iif,Vf)  must  equal  zero. 

For  developable  surfaces,  the  unit  normal  vector  N  is  constant  along  a  generator.  There- 
fore, the  rate  of  change  of  the  unit  normal  vector  in  the  u  direction  must  also  be  constant 
along  a  generator.  This  leads  us  to  the  fact  that  Nu  is  not  only  zero  at  (uf,Vf)  but  also 
zero  along  the  u  =  Uf  isoparametric  line.  Therefore,  for  a  given  u  —  Uf,  equation  (7.3) 
becomes 

L(uf,v)  =  -Ruu(w/,u)  •  N(ti/,u)  =  0  (7.6) 

for  0  <  v  <  1.  Consequently,  the  entire  generator  consists  of  a  line  of  flat  points.  I 

For  a  developable  surface,  the  inflection  line  is  a  generator  which  consists  of  a  line  of 
flat  points  and  coincides  with  a  change  in  sign  of  the  nonzero  principal  curvature.  In  other 
words,  the  nonzero  principal  curvature  has  a  different  sign  on  either  side  of  an  inflection 
line. 


A  developable  surface  can  not  possess  spherical  umbilics  since  one  of  the  principal  curvatures  is  always 
zero. 
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The  inflection  line  can  be  detected  by  finding  u  —  u^  such  that  L(v,f,vn)  =  0  where  vn 
is  an  arbitrary  constant  between  0  and  1.  L(u,vn)  =  0  can  be  written  as 

Ruu{U,Vn)  •  —— - — 77=0-  (7.7) 

IRu^Un)  X  R,,(ti,Un)| 

Since  we  are  assuming  a  regular  surface  such  that  |Ru  x  Rv|  ^  0,  we  need  only  set  the 
numerator  of  equation  (7.7)  to  zero.  Thus, 

|Ruu(«,vn)      Ru{u,vn)     R„(u,v„)|=0. 

For  a  polynomial  surface  with  degree  n  in  the  u  direction,  this  results  in  a  univariate 
polynomial  equation  of  degree  (3n  —  4)  in  u 

\DuZv        ZuVvJ-Euu        \%vZu        •^u^vlUuu  T  \%uDv        ^vVuJ^uu  =  "•  \>  •&) 

If  the  surface  is  expressed  in  a  piecewise  polynomial  form  such  as  a  B-spline  representation, 
equation  (7.8)  must  be  applied  to  each  polynomial  segment  separately.  The  univariate  poly- 
nomial equation  can  be  robustly  and  efficiently  solved  by  the  Interval  Projected  Polyhedron 
algorithm  [24,  39].  Once  we  find  the  u  value,  we  check  if  the  nonzero  principal  curvature 
changes  sign.  If  the  sign  does  not  change,  the  generator  is  not  an  inflection  line.  A  simple 
example  is  given  by  H(u,v)  =  (u,v,u*).  Since 

2^        __,        ^        (-4U3,0,1) 


G=l,         Run  =  (0,0, 12uJ)      and      N  = 


VlQu6  +  1  ' 
we  have 

-12u2 


K     = 


v^!6w6  +  1 


It  is  apparent  that  k*  becomes  zero  when  u  =  0;  however  k*  does  not  change  sign  when  u 
moves  from  u  <  0  to  u  >  0.  Therefore,  u  =  0  is  not  an  inflection  line. 

At  flat  points  the  principal  directions  are  indeterminate  and  the  orthogonal  net  of  lines 
of  curvature  may  have  singular  properties.  In  the  following  we  investigate  the  pattern  of 
the  lines  of  curvature  near  the  line  of  non-generic  flat  points.  Locally  any  surface  is  the 
graph  of  a  differentiable  function  [8].  Given  a  point  p  of  a  surface  S  we  can  choose  the 
coordinate  axis  so  that  the  origin  O  of  the  coordinate  is  p  and  the  z  axis  is  directed  along 
the  positive  normal  of  S  at  p;  thus,  the  xy  plane  coincides  with  the  tangent  plane  of  S  at 
p.  It  follows  that  a  neighborhood  of  p  in  S  can  be  represented  in  the  form  z  —  h(x,  y)  with 
/i(0,0)  =M0,0)  =M°>°)  =0. 

By  taking  into  consideration  that  h(0, 0)  =  hx(0, 0)  =  hy(0, 0)  =  0,  the  Taylor  expansion 
of  the  z  component  of  the  surface  becomes 

h{x,y)     =     ^[^2^xx(0,0)+2a;y/ixy(0,0)+y2^y(0,0)] 

+  y[x3hxxx(0,0)  +  3x2yhxxy{0,0)  +  3xy2hxyy{0,0)  +  y3/iyyy(0,0)] 
+JR(x,y)(|x,y|3)  (7.9) 


where  R(x,  y)  is  a  remainder  term  with  limx_,.o,y->o  R(x,  y)  =  0  and  |x,  y\  =  y/x2  +  y2. 
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This  local  approximation  will  now  be  applied  to  developable  surfaces. 

Lemma  7.2.1  A  developable  surface  is,  in  general,  locally  a  parabolic  cylinder  and  becomes 
a  cubic  cylinder  at  the  inflection  line. 

Proof:  Let  us  consider  an  orthogonal  Cartesian  reference  frame  attached  to  the  surface  at 
an  arbitrary  point  (uq,vo).  We  choose  unit  vectors  £*-r  x  N,  t^t  and  N  as  the  directions 
of  x,  y  and  z  axes  such  that  the  y  axis  coincides  with  the  generator  u  =  uq,  the  z  axis 
coincides  with  the  normal  vector  and  the  x  axis  is  orthogonal  to  both  axes.  Equation  (7.9) 
reduces  to  [32] 

h(x,y)  =  -x2hxx{0,0)  +  -[x3hxxx(0,0)  +3x2yhxxy{0,0)]  +  R(x,y)(\x,y\3).        (7.10) 

Using  the  inverse  function  theorem,  hxx(0,0),  hxxx(0,0)  and  hxxy(0,0)  can  be  obtained  as 
a  function  of  u  and  v.  At  a  the  line  of  inflection,  hxx(0,0)  and  hxxy(0,0)  become  zero  [32], 
and  equation  (7.10)  reduces  to 

3 


h(x,  y)  =  i(N  •  R„uu)     J EQG_  p2  \ 

where  the  coefficient  of  a;3  is  evaluated  at  (mo,^o)-  From  equation  (7.10)  it  is  apparent  that 
for  small  x,  the  quadratic  term  dominates.  From  equation  (7.11)  it  is  apparent  that  at  an 
inflection  line,  the  surface  is  locally  a  cubic  cylinder.  I 

Theorem  7.2.2  There  is  only  one  line  of  curvature  that  passes  through  each  flat  point  on 
an  inflection  line,  and  that  line  of  curvature  is  orthogonal  to  the  direction  of  the  generator. 


Proof:  By  Lemma  7.2.1  the  developable  surface  is  expressed  locally  as  a  cubic  cylinder. 
If  we  rewrite  the  cubic  term  of  equation  (7.11)  in  terms  of  polar  coordinates  by  substituting 
x  =  r  cos  9  we  obtain 

h(e)  =  j(N.ruuu)U EQG_ FA    cos39.  (7.12) 

Since  h(0  +  rr)  =  —h(6),  h(6)  is  an  antisymmetric  function  of  9.  The  roots  of  ^|  =  0  will 
give  the  angles  where  local  maxima  and  minima  of  h(9)  may  occur  around  the  flat  point. 
The  equation  can  be  restricted  to  the  range  0  <  9  <  2n  without  loss  of  generality.  The 
roots  are  easily  computed  as  9  =  0,  |,  n  and  3y.  Only  9  =  0  and  9  =  n  (which  coincide 
with  the  local  x  axis)  give  extrema,  since 

#h{i)  _^(¥)=fl 
d92  d92 

and  thus  9  =  |  and  9  =  4p  (which  coincide  with  the  local  y  axis)  provide  neither  a  maximum 
nor  a  minimum.  In  other  words,  the  inflection  line  is  not  a  line  of  curvature.  Consequently, 
there  is  only  one  line  of  curvature  (minimum  curvature  where  9  =  0  and  maximum  curvature 
where  9  =  ir)  that  passes  through  the  flat  point.  I 
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Recall  from  Section  2.3  that  the  principal  curvatures  can  be  found  by  solving  the  two 
simultaneous  equations  (2.19)  and  (2.20), 

(L  +  KE)du  +  {M  +  KF)dv    =    0 
(M  +  KF)du  +  {N  +  nG)dv    =    0, 


for  k.  Since  M  =  N  =  0,  equations  (2.19)  and  (2.20)  reduce  to 

(L  +  nE)du  +  nFdv     =    0     and 
nFdu  +  nGdv     =    0. 


(7.13) 
(7.14) 


We  can  trace  the  lines  of  curvature  which  pass  through  the  flat  points  of  an  inflection  line 
by  providing  the  starting  points  and  shifting  outwards  in  the  directions  0  and  n  from  the 
flat  points  or,  equivalently,  along  the  positive  and  negative  local  x  axis.  The  tracing  can 
be  carried  out  by  integrating  the  initial  value  problem  for  the  coupled  ordinary  differential 
equations  (7.13)  and  (7.14). 

Generic  umbilics  (including  flat  points)  are  stable  with  respect  to  small  perturbations  of 
the  function  representing  the  surface,  while  non-generic  ones  are  unstable  [32].  In  generic 
cases  umbilics  are  isolated  [20];  thus  the  inflection  line,  which  consists  of  line  of  flat  points,  is 
non-generic  and  therefore  unstable.  In  the  following  we  give  a  couple  of  numerical  examples 
that  demonstrate  the  instability  of  the  line  of  flat  points  along  the  inflection  line  with  respect 
to  perturbations. 

The  example  surface  is  a  degree  (3-1)  integral  Bezier  patch  which  is  constructed  by  the 
method  described  in  Section  5.3,  see  Figure  7-1.  The  control  points  are  shown  in  Table  7.1. 
The  surface  has  an  inflection  line  at  u  =  0.5754,  which  has  been  computed  by  solving  the 


Po,o 

(0.0000,  0.0000,  0.0000) 

Po,i 

(0.5000,  0.0000,  2.0000) 

Pi,o 

(1.8000,  3.0000,  0.0000) 

Pi,o 

(1.8950,  2.3250,  2.0000) 

P2,0 

(3.3000,-2.0000,  1.5000) 

P2,0 

(3.0575,-1.5500,  3.1625) 

P3,o 

(4.0000,  0.0000,  0.0000) 

P3,0 

(3.6000,  0.0000,  2.0000) 

Table  7.1:  Control  Points  for  a  Degree  (3,1)  Developable  Surface 

degree  5  univariate  polynomial  equation  (7.8).  The  inflection  line  is  shown  in  Figure  7-1 
with  a  dash  dotted  line.  This  surface  has  a  net  of  lines  of  curvature  which  is  shown  in 
Figure  7-2 (a).  Solid  lines  represent  the  lines  of  maximum  principal  curvature,  while  dotted 
lines  represent  the  lines  of  minimum  principal  curvature.  The  inflection  line  is  depicted  with 
a  dash  dotted  line.  Figure  7-2 (b)  shows  a  magnification  near  the  inflection  line.  We  can 
observe  that  there  is  only  one  line  of  curvature  that  passes  through  a  flat  point  orthogonal 
to  the  inflection  line. 

We  gradually  perturb  the  control  points  of  the  surface  and  observe  the  behavior  of  the 
lines  of  curvature  which  pass  through  the  inflection  line.  The  control  points  are  perturbed  in 
the  following  manner.  Since  the  example  is  a  degree  (3-1)  patch,  it  has  8  control  points.  Each 
control  point  consists  of  three  Cartesian  coordinates  x,  y,  z;  hence  there  are  24  components 
to  be  perturbed.  A  random  number  which  varies  from  —1  to  1  is  used  to  determine  the 
perturbation  of  each  component.  Let  us  denote  the  randomly  chosen  numbers  for  each 
control  point  as  (ef  ,  e^  ,  e*A  0  <  i,j  <  3.   We  normalize  the  vector  and  add  it  to  each 
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Figure  7-1:  Developable  surface  and  its  control  polyhedron  with  inflection  line. 


inflection  line 


Figure  7-2:  (a)  Lines  of  curvature  of  developable  surface  with  inflection,  (b)  Magnification 
near  inflection  line. 
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Figure  7-3:   (a)  Lines  of  curvature  on  perturbed  surface  £  =  0.02.   (b)  Magnification  near 
u=0.57. 


Figure  7-4:    (a)  Lines  of  curvature  on  perturbed  surface  ( 
u=0.57. 


=  0.08  (b)  Magnification  near 
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control  point  as  follows: 


,'ty ,  <7-i5> 

~.y  *■  i  „z  i 


v/eu2  +  4" 


where  C,  =  0.02.  We  gradually  increase  the  perturbation  by  increasing  £  from  0.02  to  0.08 
in  steps  of  0.02. 

Figures  7-3  and  7-4  illustrate  the  behavior  of  the  lines  of  curvature  when  the  control 
points  are  perturbed.  We  can  see  from  the  figures  that  the  entire  inflection  line,  which 
consists  of  a  line  of  flat  points,  disappears.  Hence  there  is  no  singularity  in  the  net  of  lines 
of  curvature  when  a  perturbation  is  induced.  The  nonzero  principal  curvatures  on  both 
sides  of  the  former  inflection  line  2  meet  at  right  angles  near  the  former  inflection  line  and 
make  a  very  sharp  change  in  direction  (almost  a  right  angle). 


2 Once  the  control  points  are  perturbed  both  principal  curvatures  may  not  be  nonzero,  but  here  we  are 
referring  to  the  nonzero  principal  curvature  before  perturbation. 


Chapter  8 

Engineering  Example 


The  final  example  is  the  design  of  a  small  boat  using  B-spline  developable  surfaces.  The 
boat  was  designed  with  a  hard  chine  which  separates  the  bottom  and  side  of  the  boat  and 
demarks  a  tangent  discontinuity  in  the  hull  of  the  boat.  As  an  initial  stage  of  the  design,  the 
following  three  curves,  which  define  the  desired  characteristics  of  the  hull,  were  determined 
using  a  model  of  a  planing  boat: 

•  sheer  curve  -  the  line  that  runs  along  the  deck  edge  from  bow  to  stern 

•  base  curve  -  the  line  that  runs  along  the  keel  of  the  boat  (centerline  from  bow  to 
stern) 

•  chine  curve  -  the  line  that  runs  from  the  bow  to  the  lower  corner  of  the  stern 

Figure  8-1  shows  the  terminology  used  in  this  section.  The  coordinate  system  is  located 
so  that  the  origin  is  located  on  the  baseline  at  the  bow  of  the  boat  with  the  x  axis  pointing 
aft,  the  y  axis  pointing  to  starboard  and  the  z  axis  pointing  up.  (Note  that  the  baseline, 
a  zero  reference  line,  is  different  from  the  base  curve  defined  above.)  Due  to  the  symmetry 
of  the  boat,  we  design  only  the  starboard  side  and  reflect  those  values  across  the  plane  of 
symmetry  to  obtain  the  port  side.  The  cubic  B-spline  representation  of  the  three  curves  is 
shown  in  Table  8.1. 


Transom 


Chine  Sheer 


Figure  8-1:  Planing  Boat  Terminology 


Design  Using  3D  Directrices 

Employing  the  design  philosophy  described  in  Chapter  5,  the  chine  curve  was  designated 
as  the  design  curve  and  the  sheer  and  base  curves  as  target  curves,  i.e.  r#(u),  to  construct 
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Sheer  Curve 

Chine  Curve 

Base  Curve 

Control 
Points 

(0.00,  0.00,  9.00) 
(6.86,  7.10,  8.22) 
(21.6,  8.93,  6.25) 
(36.9,  8.73,  5.86) 
(45.0,  7.65,  6.10) 

(1.40,  0.00,  5.30) 
(10.5,  7.53,  1.93) 
(25.7,  7.85,  1.28) 
(40.4,  7.46,  1.27) 
(44.1,  7.20,  1.70) 

(1.40,  0.00,  5.30) 
(2.26,  0.00,  -0.21) 
(22.6,  0.00,  -0.10) 
(36.3,  0.00,  -0.10) 
(44.1,  0.00,  0.50) 

Knot 
Vectors 

0,  0,  0,  0,  .5, 
1,  1,  1,  1 

0,  0,  0,  0,  .75, 
1,  1,  1,  1 

0,  0,  0,  0,  .5, 
1,  1,  1,  1 

Table  8.1:  Target  Curve  Data 


the  side  and  bottom  respectively.  Initially,  the  endpoints  of  the  sheer  and  base  curves  were 
used  to  construct  the  side  and  bottom  of  a  single  surface  each;  however,  this  resulted  in  a 
flat  side  and  bottom.  To  remedy  this  situation,  the  chine  was  split  into  two  B-spline  curves 
at  a  parameter  value  of  u  =  0.5,  and  one  additional  knot  was  inserted  at  the  midpoint  of 
the  forward  B-spline  curve.  The  resulting  curves  are  shown  in  Table  8.2.  Thus,  the  side 
and  bottom  are  constructed  of  two  surfaces  each,  allowing  more  flexibility  in  the  design, 
but  losing  the  C2  continuity  along  the  surface.  The  surfaces  are,  however,  at  least  Gl 
continuous  between  the  forward  and  aft  portions.  Since  the  design  curve  was  split  into  two 
curves  that  maintain  C2  continuity  between  them,  the  tangent  at  the  connecting  point  is 
continuous.  Since  both  patches  are  developable,  both  surfaces  share  the  same  tangent  plane 
at  the  connecting  generator.  Hence,  the  patches  are  Gl  continuous. 

Note  that  the  forward  and  aft  surfaces  join  along  a  generator  which  is  a  straight  line. 
The  planar  transom  (stern)  joins  the  aft  surfaces  along  generators  as  well. 

The  endpoints  for  the  side  and  bottom  surfaces  were  selected  so  as  to  force  the  resulting 
directrix  as  close  to  the  target  curves  as  possible  while  maintaining  a  maximum  Gaussian 
curvature  on  the  order  of  10-5.  The  data  for  each  developable  surface  on  the  starboard  side 
of  the  boat  is  listed  in  Table  8.2.  Measurements  are  in  feet.  The  port  side  of  the  boat  has 
the  same  measurements,  with  a  negative  y  value. 


Patch 

Design  Curve 
Control  Points 

Knot 
Vector 

Endpoints 

Plane 
Definition  Point 

Gaussian 
Curvature 

Forward 
Side 

(1.40,  0.00,  5.30) 
(4.17,  2.32,  4.26) 
(10.7,  5.86,  2.56) 
(18.2,  7.09,  1.87) 
(21.9,  7.38,  1.66) 

0,  0,  0,  0 

0.5 
1,1,  1,  1 

(0.00,  0.00,  9.00) 
(20.0,  8.30,  6.74) 

(20.0,  0.00,  6.74) 

8.18  x  10~5 

Forward 
Bottom 

(1.40,  0.00,  5.30) 
(14.0,  0.00,  0.00) 

(0.00,  0.00,  0.00) 

0.00 

Aft 
Side 

(21.9,  7.38,  1.66) 
(25.6,  7.67,  1.44) 
(33.0,  7.66,  1.28) 
(40.4,  7.46,  1.27) 
(44.1,  7.20,  1.70) 

0,  0,  0,  0 
0.5 

1,  1,  1,  1 

(20.0,  8.30,  6.74) 
(45.0,  7.70,  6.10) 

(45.0,  0.00,  6.10) 

6.62  x  10~5 

Aft 
Bottom 

(14.0,  0.00,  0.00) 
(44.1,  0.00,  0.50) 

(0.00,  0.00,  0.00) 

9.74  x  10_b 

Table  8.2:  Developable  Surface  Data,  Starboard  Side 
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Constrained  by  Plane 

All  the  second  directrices  were  constrained  so  that  they  lie  in  a  plane  defined  by  the  two 
designated  endpoints  of  rs(u)  and  a  third  point  g  provided  by  the  user.  This  point  g  is 
termed  the  Plane  Definition  Point  and  is  shown  in  Table  8.2.  For  the  bottom  surfaces, 
the  constraining  plane  was  the  xz-plane  which  divides  the  boat  into  symmetric  halves.  For 
the  side  surfaces,  the  constraining  plane  was  positioned  so  that  a  line  in  the  plane  and 
perpendicular  to  the  x  axis  would  be  parallel  to  the  xy  plane. 

Figure  8-2  shows  the  resulting  developable  surfaces  with  the  design  lines  shown  as  dotted 
lines  for  reference.  Note  that  the  forward  portion  of  the  bottom  of  the  boat  is  a  triangular 
degenerate  patch.  Figure  8-3  is  a  bottom  view  of  the  boat.  One  side  shows  the  rulings  of 
the  developable  surfaces  while  the  other  side  is  a  shaded  view. 


Figure  8-2:  Boat  Composed  of  Developable  Degree  (3-1)  B-Spline  Surfaces 


Figure  8-3:  Boat  Bottom  View 
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Development 

The  plane  development  of  each  surface  of  the  boat  was  computed,  thus  providing  cutting 
information  for  the  flat  steel.  Using  the  method  described  in  Section  5.4,  the  planar  surfaces 
shown  in  Figure  8-4  were  constructed.  All  developments  were  within  a  tolerance  of  10~6. 


(a)  Aft  Side 


(b)  Forward  Side 


(c)  Aft  Bottom 


(d)  Forward  Bottom 


Figure  8-4:  Boat  Surfaces  Developed  onto  a  Plane 


Geodesies 


At  the  Electric  Boat  Corporation  shipyard,  many  developable  surfaces  are  used  in  the 
construction  of  submarines,  since  a  large  portion  of  the  hull  is  cylindrical.  When  bending 
the  steel  into  shape,  accuracy  is  checked  with  detailed  measurements  of  geodesies  drawn  on 
the  surface.  When  the  surface  is  flat,  straight  lines  are  drawn  between  the  corners  of  the 
surface,  forming  an  X  as  shown  in  Figure  8-5(a).  The  position  of  this  line  on  the  curved 
surface  is  easily  computed  using  the  method  described  in  Chapter  6.  As  an  example,  the 
geodesies  on  the  Forward  Side  section  have  been  calculated.  The  uv  parameter  values  are 
shown  in  Figure  8-5(b)  and  the  geodesies  on  the  developable  surface  are  shown  in  Figure  8- 
5(c). 

Lines  of  Curvature 

These  surfaces  have  no  inflection  lines,  but  the  lines  of  curvature  were  calculated  as  described 
in  Chapter  7  and  are  shown  in  Figure  8-6.  The  dashed  lines  represent  lines  of  minimum 
curvature  and  the  solid  lines  are  lines  of  maximum  curvature.  Note  that  one  set  is  always 
composed  of  straight  lines  that  correspond  to  the  generators. 
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(a)  Surface  Developed  onto  Plane 


(b)  uv  Parametric  Space 


(c)  Geodesic  on  3D  Surface 


Figure  8-5:  Geodesies  on  the  Forward  Side  Section 
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(b)  Forward  Side 


(a)  Aft  Side 


(d)  Forward  Bottom 


(c)  Aft  Bottom 


Figure  8-6:    Lines  of  Curvature  on  Boat  Surfaces  (Solid  Lines  are  Maximum  Curvature 
Lines,  Dashed  Lines  are  Minimum  Curvature  Lines) 
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Chapter  9 

Conclusions  and  Recommendations 


9.1  Conclusions 

In  this  thesis  a  novel  method  to  design  complex  objects  in  terms  of  B-spline  developable 
surfaces  focusing  on  user  friendliness  is  developed.  The  user  only  needs  to  specify  the 
design  B-spline  curve  and  the  two  end  rulings,  and  the  developable  surface  in  a  B-spline 
representation  is  automatically  generated.  Although  optimization  techniques  are  relied 
upon,  in  most  cases  the  computed  results  have  a  Gaussian  curvature  on  the  order  of  10 
or  less  unless  the  design  of  a  very  complex  surface  is  attempted.  The  computational  time 
for  the  design  of  all  the  examples  in  this  thesis  is  on  the  order  of  a  few  tenths  of  seconds  on 
a  100  MHz  Silicon  Graphics  workstation. 

In  addition,  two  differential  geometry  properties  of  developable  surfaces  which  have 
practical  engineering  applications  were  investigated.  The  first  property  was  the  definition 
of  inflection  lines  on  developable  surfaces,  which  are  critical  in  the  manufacture  of  such 
surfaces.  We  proved  that  an  inflection  line  on  a  developable  surface  occurs  on  a  generator 
which  consists  of  entirely  flat  points,  and  derived  a  single  variable  nonlinear  polynomial 
equation  to  locate  it.  It  was  also  found  that  for  each  flat  point  on  an  inflection  line  of  a 
developable  surface,  there  is  only  one  line  of  curvature  that  passes  through  it,  and  that  line 
of  curvature  is  orthogonal  to  the  direction  of  the  generator. 

The  second  differential  geometry  property  investigated  was  the  location  of  geodesies 
between  two  points  on  a  developable  surface.  It  was  shown  that  the  problem  can  always  be 
converted  from  a  boundary  value  problem  to  an  initial  value  problem  with  a  corresponding 
savings  in  time  for  the  solution. 

Finally,  these  techniques  were  applied  to  a  boat  model,  a  complex  object  that  was 
represented  entirely  with  B-spline  approximate  developable  surfaces  with  cubic  directrices. 
The  Gaussian  curvature  was  on  the  order  of  10-5.  Geodesies  and  lines  of  curvature  were 
calculated  on  these  surfaces  as  well. 

9.2  Future  Work 

The  ultimate  goal  of  this  research  is  to  be  a  useful  tool  in  the  design  and  manufacturing  of 
ship  hulls.  Therefore,  future  research  should  take  this  work  in  that  direction.  Some  specific 
areas  of  possible  research  are  described  below. 

Approximation.  The  approximation  of  free  form  surfaces  by  developable  surfaces  would 
be  extremely  helpful  in  converting  existing  ship  designs  to  ones  including  more  developable 
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surfaces.  Some  preliminary  research  has  been  completed  by  Elber  and  Cohen  [9]  into  repre- 
senting surfaces  as  ruled  surfaces,  which  could  possibly  be  extended  to  include  developable 
surfaces.  In  addition,  research  on  approximation  using  developable  surfaces  was  conducted 
by  Hoschek  and  Pottmann  [23]  using  the  representation  of  developable  surfaces  as  the  en- 
velope of  a  family  of  planes.  Alternatively,  Cho  [6]  discussed  the  approximate  development 
of  non-developable  doubly  curved  surfaces  onto  a  plane,  which  may  also  be  useful  in  the 
above  approximation  problem. 

Manufacturing.  Once  the  surface  is  developed  into  a  plane  and  the  steel  is  cut  to  the 
correct  size,  it  must  be  bent  into  the  proper  shape.  Additional  research  should  investigate 
the  proper  amount  of  curvature  to  be  induced  by  rolling  or  pressing  and  instructions  for 
the  shop  floor  should  be  output.  Although  the  lines  of  curvature  are  instrumental  in  this 
bending,  mechanical  considerations  such  as  springback  must  be  included  in  the  design  as 
well.  Much  research  has  been  done  in  this  area  by  Hardt  et  al  [19,  18,  17].  The  ultimate 
goal  is  for  the  engineer  to  design  the  ship  hull  and  output  information  that  is  directly  useful 
to  the  workers  on  the  shop  floor.  To  this  end,  the  operating  characteristics  of  bending 
equipment  should  be  studied  and  properly  accounted  for. 

Geometry.  Several  areas  of  interest  remain  for  investigation  within  the  differential  ge- 
ometry discussed  in  this  thesis.  First,  a  suitable  non-dimensional  Gaussian  curvature  ex- 
pression should  be  developed  to  better  judge  the  developability  of  surfaces  designed  using 
the  minimization  process  described  in  Section  5.3. 

Although  strip  surfaces  of  any  degree  can  be  handled  by  the  computer  programs  devel- 
oped during  this  research,  the  programs  for  the  surfaces  with  arbitrary  3D  directrices  can 
handle  only  B-spline  developable  surfaces  with  cubic  directrices.  It  may  be  beneficial  to 
investigate  higher  degree  surfaces  which  may  include  more  flexibility  and  would  be  able  to 
more  closely  approximate  the  desired  output  surface. 

In  addition,  the  surfaces  in  this  thesis  are  all  integral  B-spline  surfaces.  Extension 
to  include  rational  B-spline  developable  surfaces  such  as  those  investigated  by  Pottmann 
and  Farin  [37]  and  Lang  and  Roschel  [29]  would  allow  the  representation  of  more  complex 
surfaces. 

One  difficulty  in  designing  surfaces  using  this  method  is  that,  for  surfaces  with  3d 
directrices,  there  is  no  method  to  determine  in  advance  if  the  surface  will  include  the  edge 
of  regression.  Therefore,  regular  surfaces  are  achieved  through  trial  and  error,  and  hence 
regularity  conditions  should  be  investigated  more  fully. 

Design.  Much  research  has  been  conducted  in  design  for  producibility  in  ship  hulls.  For 
example,  Lamb  [27]  shows  a  method  to  construct  a  bulbous  bow  from  regular  surfaces  and 
discusses  the  simplification  of  bow  and  stern  sections  to  include  more  developable  surfaces. 
This  should  be  expanded  to  investigate  the  design  of  an  entire  ship  hull  from  developable 
surfaces.  Once  a  ship  hull  is  designed,  it  should  be  hydrodynamically  tested  to  compare  the 
performance  to  a  similar  hull  that  is  not  constructed  of  developable  surfaces.  The  relative 
magnitude  of  the  savings  due  to  producibility  and  the  losses  due  to  performance  should  be 
compared. 
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